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ABSTRACT 

We present a detailed study of the analytic structure, BPS spectra and superconformal 
points of the N — 2 susy SU(2) gauge theories with Nf = 1, 2, 3 massive quark hypermultiplets. 
We compute the curves of marginal stability with the help of the explicit solutions for the low 
energy effective actions in terms of standard elliptic functions. We show that only a few of 
these curves are relevant. As a generic example, the case of Nf = 2 with two equal bare 
masses is studied in depth. We determine the precise existence domains for each BPS state, 
and show how they are compatible with the RG flows. At the superconformal point, where two 
singularities coincide, we prove that (for Nj = 2) the massless spectrum consists of four dist inct 
BPS states and is 5-invariant. This is due to the monodromy around the superconformal point 
being S, providing strong evidence for exact 5-duality of the SOFT. For all Nf, we compute 
the slopes u of the /3-functions at the fixed point couplings and show that they are related to 
the anomalous dimensions a of u = (tr</> 2 ) by u — 2(a — 1). 


* Unite Propre de Recherche 701 du CNRS, associee a l’Ecole Normale Superieure et a l’Universite 
Paris-Sud. 



1. Introduction 


Many new insights into the physics of strongly coupled gauge theories have been obtained 
by the study of N = 2 supersymmetric Yang-Mills theories [1,2]. Two particularly interesting 
phenomena that occur in these four-dimensional theories are the discontinuities of the BPS 
spectra [3,4,5] on the Coulomb branch and the occurrence of superconformal points that are 
believed to lead to non-trivial interacting 4D superconformal held theories (SOFT) [6,7]. 

In the present paper we study the probably simplest case where both phenomena occur: 
the SU(2) theories with massive quark hypermultiplets. The discontinuities of the BPS spectra 
and the properties at the superconformal points are of course closely related and our study of 
the first will shed important new light on the second, and vice versa. 

The main ingredients in the study of the N — 2 gauge theories are duality and holomor- 
phicity of the low-energy effective action. The properties of the latter are encoded in a (hyper) 
elliptic curve which also determines the abelian charges of the theory. The BPS states consti¬ 
tute a particularly important sector of the Hilbert space: all perturbative and known solitonic 
states are BPS states. Their masses are determined by the abelian charges and hence by the 
(hyper) elliptic curve. Quite surprisingly, the latter also allows us to extract and determine the 
existence domains of the BPS states on the Coulomb branch of the moduli space A4. Indeed, 
BPS states generically are stable, except on real codimension one hypersurfaces in Ni (i.e. real 
dimension one curves for SU(2)). These instability hypersurfaces are determined in terms of 
the abelian charges and the quantum numbers of the BPS state. 

For the SU(2) theories without or with Nf massless quark hypermultiplets (Nf = 1,2,3) 
there only is a single instability curve on M ~ C. This curve is closed and goes through 
the singular points on M.. The BPS spectra are different inside and outside the curve [3,4,5]. 
We found that in the region outside the curve, part of which is the semiclassical domain, all 
semiclassically expected states exist. These are: the dyons (n e ,n m ) = (2n, 1) for Nf = 0, 
the quarks (1,0) and the dyons (n, 1) for Nf = 1,2,3, with in addition dyons (2 n + 1,2) of 
magnetic charge two for Nf = 3. There also is the W-boson (2,0) for all Nf. In the region 
inside the instability curve, which always is a region of strong coupling, there only exist the 
states that are responsible for the singularities, i.e. the states that can become massless. For 
Nf = 0 e.g., these are the dyons (0,1) and (—2e, 1) while for Nf = 2 these are (0,1) and 
(—e, 1). Here e equals +1 in the upper and —1 in the lower half plane. Almost all BPS states 
disappear when crossing the curve from outside to inside: they “decay” into the only two (or 
three) existing states. In particular, there is no W-boson inside the curve, although the gauge 
symmetry is always broken SU(2) —» U(l). 

The determination of these BPS spectra was relatively easy since there was only one 
instability curve. The generic situation of SU(n) is much more complicated: each BPS state 
has its own family of possible instability hypersurfaces. At first sight it seems hopeless to study 
this general case by the methods of [3,4], However, a similar situation already occurs in the 
SU(2) theories with massive hypermultiplets. This is the case we study in the present paper. It 
will turn out that among the multitude of possible decay curves only a relatively small subset 
is relevant, and one obtains a very clear picture of the existence domains of the various BPS 
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states. We are quite confident that this structure can be exploited to also determine the exact 
BPS spectra at any point on the Coulomb branch of the SU(n) moduli space. 

Another feature, absent in the massless SU(2) theories, but present for the higher SU(n) 
theories or the SU(2) theories with massive hypermultiplets, is the existence of snperconformal 
points. These can occur at special points on the Coulomb branch in the SU(n) theories where 
singular lines intersect, or in the massive SU(2) theories at special values of the masses when 
singular points collide. In any case, at a superconformal point two or more mutually non-local 
BPS states simultaneously become massless [6]. To study these points, it will turn out to be 
most fruitful to view them as resulting from such a coincidence of individual singularities. 

In this paper, we will discuss the BPS spectra of the SU(2) theories with Nf = 1,2,3 
massive quark hypermultiplets, and present an in depth study of the Nf = 2 theory with equal 
bare masses. This latter case is sufficiently generic to exhibit all interesting new phenomena. 
We determine the precise existence domains for each BPS state and confirm these results by 
many additional consistency checks. We will also see how the whole set of decay curves and 
BPS spectra very consistently behaves under the RG flow from one Nf theory to another. A 
subtlety present in the massive theories is related to the abelian .^-charges (i — 1,... Nf) [8] 
which we need to determine for all BPS states. Some indications on the BPS spectra of the 
massive theories were already obtained in [9] within the geodesics approach from string theory 
[10]. However, it only provided some partial and pointwise information on M.. 

The analytic structure and the corresponding monodromies are fundamental for our study. 
As we follow the various RG flows, we reach various superconformal points where the mon¬ 
odromies are given by the products of the individual monodromies of the coinciding singu¬ 
larities. Such a composite monodromy M sc is quite special, acting e.g. as ^-duality relating 
mutually non-local states as monopoles and quarks. We argue that it should be an exact quan¬ 
tum symmetry of the massless BPS sector and thus of the superconformal field theories. This 
monodromy (completely) characterises the SCFT allowing us to compute scaling dimensions. 
We show that the SCFT is determined in terms of a single integer k — 1, 2 or 3 characteristic 
of M sc . This is reminiscent to the study of the relevant deformations of the singular curves 
y 2 = x 3 in [7]. As we follow the RG flow further, the singularities separate again but the 
analytic structure is changed, providing us with an explanation of how the nature of certain 
singularities can change from being due to a massless quark at weak coupling to being due to 
a massless dyon at strong coupling. 

Let us outline the organisation of the present paper and some of our results: first, in Section 
2, we recall the elliptic curves for the massive Nf = 1,2,3 theories and discuss the quantum 
numbers of the BPS states associated with the singularities and how they change under the 
various RG flows one can study. This brings us naturally to a discussion of the superconformal 
points and how they are classified by an integer k — 1, 2, 3 through their monodromies. Then, 


v Where comparable, both results agree. Although there is a slight discrepancy with the published version 
of [9], after a first circulation of the present paper, the authors of [9] have informed us that this is only 
due to some error when writing up their paper but that they actually agree with our results, see Section 
4 below. 
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we compute the basic functions cld{u ) and a(u), which are the period integrals of a certain 
meromorphic one-form, in terms of standard elliptic integrals in a form immediately suited for 
numerical computations. Technical details are postponed to the appendices. Similar results 
were also obtained independently and published recently in [11,12], In Section 3, we discuss 
general features of the BPS spectra and the decay curves for all the massive Nf = 1,2,3 
theories. In particular, we discuss the maximal possible set of BPS states compatible with the 
RG flow. 

Section 4 then is an in depth study of the Nf = 2 theory with two equal bare masses 
m. There exists a superconformal point when m — 4p, and we discuss separately the cases 
rn < 4^ and rn > 4p. We find that all curves and existence domains of BPS states have 
perfectly smooth RG flows. In particular, the crossover from small mass (m < 4jr) to large 
mass (m > 4jr) is perfectly smooth, except, in a certain sense, at the superconformal point w* 

itself. Also, as m becomes very large, in a basis op, a that flows to the a (Ab=0) only 

the states with s = 0 survive under the RG flow m —> oo to the pure gauge theory Nf = 0. 
At a given fixed point u G A4, almost all s / 0 dyons disappear, already at finite m, because 
they are “hit” by their corresponding decay curves that move outwards (to large \u\) as m is 
increased. Among the 0 states only the quarks and some special dyons exist for all finite 
m, but as rn —> oo, their BPS masses diverge and thus they simply drop out of the spectrum 
for this reason. On the other hand, the different s = 0 states decay on one and the same curve 
which flows to the curve of the Nf = 0 theory. Thus we are able to see the flow to the Nf = 0 
spectra in full detail. For each case, rri < 4^ and m > 4jr, we first discuss the general picture, 
which is then established by considering each class of BPS states separately. 

The BPS spectra for Nf — 2 , m — 4^ are then simply obtained as the limit of those for 
either rn > 4^ or rn < 4r• Nothing changes dramatically for rn — 4p, except at the value of 
u equal to the two coinciding singularities which is the superconformal point. In Section 5, 
we discuss the physics of this superconformal point in some detail. There, quite remarkably, 
one has four massless states, namely the quark (n e , n m ) s = (1,0)i and the monopole (0, l)o 
which are flavour doublets, and the two dyons (1, l)i and (—1,1)—i which are flavour singlets. 
This spectrum has an S'-duality invariance which is realized by the monodromy matrix M sc 
around the superconformal point. Since this massless sector constitutes the superconformal 
held theory, the latter should have a quantum S-duality invariance. We then discuss in general 
the Argyres-Douglas ansatz [6] for the /3-function of a theory based on the massless states at 
the superconformal point. These authors conjecture that, although there is no local action 
at our disposal to simultaneously describe all massless states, the beta might nevertheless be 
obtained by simply computing the contribution of each massless particle in the formulation of 
the theory which describes it locally and then adding together the suitably duality transformed 
contributions. We show that despite its appealing features, this ansatz is not correct, as already 
suspected by these authors because it let to irrational values for the slope u of the /3-function 
at the fixed point. We argue that the slopes oo can be computed from the low-energy effective 
actions alone and find that they are a rational numbers related to the scaling dimensions a 
of (tr0 2 ) as uo = 2(a — 1 ) . We show that this relation is in perfect agreement with N = 2 
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superconformal invariance. 

Then follow four appendices. In appendix A, we discuss the positions of the singularities 
and their different RG flows for Nj = 1,2 and 3. In appendix B, we give details on the elliptic 
integrals needed in Section 2 - and heavily used for the numerical computations of Section 4. 
In appendix C, we express the period integrals in terms of the three elliptic integrals also for 
Nf = 1,3 and for Nf = 2 with mj = m, m 2 = 0, and check the RG flows on these expressions. 
Finally, in appendix D, we study the RG flow of the Nj = 2 integrals with equal bare masses, 
thus providing some additional consistency checks. 


2. RG flows, analytic structure, 
superconformal points and period integrals 

The structure of the Coulomb branch of the asymptotically free SU(2) theories under study 
was derived in [1, 2], It is given in terms of an elliptic curve of the form 

y 2 = x 2 (x — u) + PN f (x,u,mj, A Nf ) (2.1) 


where* u = (tr0 2 ) is the gauge invariant moduli, (f) the scalar component of the N = 2 vector 
multiplet, rrij (1 < j < Nf) the bare masses of the hypermultiplets and A N f the dynamically 
generated scale of the theory. The polynomials Pj\r f are given by 


A 4 

P °=T X ’ 

p Af A? 

Pl= T m ix-—, 

P2 = - (rf (1 - u) + '-fmpnrpr - Amy + m|) , 


64 

Ao 


4 

,2 A 2 


64 


P, = -f A { x - u Y-^{x-u){ml + ml + ml) 

+ ^ (m\m\ + m\m\ + m 2 m 2 ) . 


( 2 . 2 ) 


The masses of the so called BPS states, which come in short representations of the supersym¬ 
metry algebra, can be computed using the fundamental formula 


A/bps(«) 


a/2 


n m aD{u) - n e a(u) + ^ . 


m 

! a/2 


(2.3) 


In this expression, n e and n m are two integers representing the electric and magnetic charges 
of the state, and the Sj are integers or half-integers which correspond to constant parts of the 


-k This relation is only valid up to a constant shift in the general case [13]. 
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physical baryonic charges [2,8]. If A is a meromorphic differential on the curve ( 2 . 1 ) such that 


d\ y/2 dx 
du 87 t y 

then the variables a and a 75 are given by the contour integrals 



(2.4) 


(2.5) 


for a certain homology basis ( 71 , 72 )- The BPS mass formula (2.3) is compatible with duality 
transformations of the form 


/ °D \ 


/ (, D \ 

( ^ \ 


( Ue \ 

a 

—> M 

a 

i I m 

—> M* 

IT'in 

\m/y/2) 


\ mj y/2 ) 

V a J 


V S / 


where 



/ a 

P 

A 


/ a 

p 

°\ 

M = 

7 

5 

9 

, M* = 

7 
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^0 

0 

l) 


V ag - 7 / 

Pg-Sf 1/ 


( 2 . 6 ) 


(2.7) 


with a , /3, 7 , 5 integers such that det M — 1 and /, g integers or half-integers. Note the useful 
relations 


(MiM 2 y = AffAfJ , (AT) 1 = (AT 1 )* . 


( 2 . 8 ) 


One particularly important class of duality transformations corresponds to the monodromy 
transformations a and ajj undergo when encircling a singularity in the u plane. These sin¬ 
gularities are due to dyons (n e ,n m )s d lying in a d-dimensiona.l representation of the flavour 
group becoming massless and occur when the curve (2.1) has a vanishing cycle, i.e. when 
the discriminant of the corresponding polynomial (2.2) vanishes. When all the non-zero bare 
masses are equal, which is the only case we will study in the following for the seek of simplicity, 
the mass term in (2.3) is of the form ms/V 2, and the monodromy matrices then read 



A 

dn e n m 

dn g 

—dn e s 

M (ne,nmV d 


-dn 2 m 

1 + dn m n e 

—dn m s 


V 

0 

0 


1 


A 

dn e n m 

dn g 

°\ 


M* — 

(■ n e ,nm)f d 


-dn 2 m 

1 + dn m n e 

0 



V 

—dsn m 

dsn e 

J 



When two singularities coincide, the monodromy is given by the product of two matrices of 
the type (2.9). Note that the fact that a and an pick up constants multiple of m/2y / 2 under 
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monodromy transformations is possible because the differential A has poles with residues such 
that 

2ni resA = ■, ti G -Z . (2-10) 

i v2 2 

2.1. RG flows and superconformal points 

The theories for different Nf are related by the renormalization group flow when some 
of the bare masses of the hypermultiplets are sent to infinity and the corresponding states 
can be integrated out (see Appendices A and C). For instance, if we let m 3 —> 00 while 
keeping A^ = fixed, the polynomial P 3 in (2.2) flows to P 2 ■ Another flow which we 

will study in great details in the following is mi = m 2 = m —> 00 with mA 2 = Aq fixed 
which allows to obtain the Nf = 0 theory directly from the Nf = 2 theory. A particularly 
important phenomenon that must occur during the RG flow is that the quantum numbers of 
some of the particles becoming massless change. Typically, when m —> 00 some singularities 
are at weak coupling and correspond to elementary quarks becoming massless while when 
m —> 0 these singularities move towards the strongly coupled region where only magnetically 
charged states can become massless. The transformation of the quantum numbers must be 
implemented by a SL(2, Z) matrix U. Since the electric and magnetic quantum numbers of 
the singularities at m — 0 can in general be deduced on physical grounds, by using for instance 
the discrete symmetry acting on the Coulomb branch when it exists, the general form of U 
can be determined a priori. Two cases can then arise: either the matrix U does or does not 
correspond to any monodromy matrix of the form (2.9). If it does not, we will see that (for 
real m) the quantum numbers change when at least two singularities coincide. Such points 
were discovered by inspection in [7]. The low energy theory at these points is believed to 
correspond to an interacting N = 2 superconformal theory [6,7]. We will discuss the physics 
of these theories in Section 5, but let us right now display the monodromy matrix around such 
a superconformal point, discuss some general constraints that it must satisfy in order to be 
compatible with conformal invariance, and determine the way the quantum numbers of the 
singularities change. 


• —> >n § 

M 1 M 2 M 1 M 2 

Fig. 1 : After the collision, the quantum numbers at the singularities and thus the 
associated monodromy matrices change according to equations (2.15) or (2.16). 


The generic case of two colliding singularities <j\ and <72 is represented in Fig. 1. The sin¬ 
gularity (7j whose monodromy is Mj is due to dj hypermultiplets ( Pj,qj) Sj becoming massless. 
Note that we choose a 1 to be below the cut produced by 02 - If the two cuts actually coincide 
on the left of cr 1 , this will mean that the quantum numbers (pi,qi) Sl at &l are computed by 
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looking at the solution (ajj,a) below the cut. With this convention, the monodromy at the 
superconformal point is 

M sc = M 2 M 1 . (2.11) 

The eigenvalues of this matrix depend on one symplectic invariant integer parameter 

k = did 2 (piq 2 - P 2 qi) 2 (2-12) 

which is such that tr M sc = 2 — k. We must have k > 0 since otherwise the particles becoming 
massless are mutually local and the low energy theory is simply N = 2 super QED with a 
dual photon. Superconformal invariance implies that the dimensions of a and ajj must both 
be equal to one, and thus that near the superconformal point w* the following expansion is 
valid 

a D {u) =a D (u *) + c D (u - «*) 1/a + o((u - u*) l/a ) 
a(u) —a(u *) + c(u - w*) 1 /® + o((w - w*) 1 /®), 

where cjj and c are constants and a is the anomalous dimension of the operator u. This is 
possible only if the eigenvalues of M sc are of modulus one, which rules out the cases k > 4. 
Moreover, when k — 4, M sc is, up to a global sign, conjugate to a certain non-zero power of 
T. This would be the signal of logarithmic terms in the asymptotic expansion of a and ajj, 
which again are ruled out by conformal invariance. Thus we conclude that 1 < k < 3. The 
three allowed values of k correspond to three inequivalent SOFT. We will bring some evidence 
in Section 4 that the SOFT is actually (fully) characterized by this integer k, by computing 
critical exponents as functions of k. 

Denoting with a prime the quantities corresponding to the new quantum numbers after 
the collision of the singularities, as indicated in Fig. 1, one must have 

M 2 Mi = M' 2 M[. (2.14) 

The multiplicity dj of a given singularity cry, which can be directly read off the discriminant 
of the relevant polynomial (2.1), (2.2), cannot change at the superconformal point. Thus one 
must either have dj = dj, in which case (2.14) can be solved non trivially by 

M' = UMjU~\ with U = M” 1 , (2.15) 

or d' x = o? 2 , d' 2 = d\. in which case (2.14) can be solved by 

M[ = UM 2 U~ x , M' 2 = UMiU _1 , with U = Aff 1 or U = M 2 . (2.16) 

The quantum numbers are then changed according to the matrix U*, see (2.7). When k = 3, 
there is another consistent solution to (2.14), which corresponds to U = I + M” 1 if dj = d'j or 
to U = M^ 1 + M 2 if d'i = o?2- The matrix U is then not of the form (2.7), but its 2 x 2 upper 
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left block part is still in SL(2, Z)* The corresponding matrix U* can then be straightforwardly 
deduced by equating the central charges of the original and image states. This latter case, 
with d\ = cfe, is actually realized in the Nj = 3 theory with three equal masses which we study 
in the next subsection. 


2.2. RG flows and analytic structure 


In Fig. 2, we have reprensented three different RG flows which illustrate the discussion of 
the previous subsection. 

The first case corresponds to the flow from the Nf = 2 theory to the Nj- — 1 theory obtained 
by sending m2 to infinity while keeping m2 A| = A| fixed, and m\ = 0. No superconformal 
point is needed in this case. The dyon singularity (1,1) 1 changes to one component of the 
elementary quark ( 1 , 0 ) 1 when crossing the cut produced by the ( 0 , l)o singularity as m 2 is 
increased. In order to recover a more conventional analytic structure in the m 2 —> 00 limit, one 
can move the cut originating at (1, 0)i to the right as indicated in the Figure. Physically, this 
simply amounts to shifting the 6 angle by an integer multiple of 2n in the upper half u-plane, 
which is an unphysical transformation. Mathematically, one keeps the solution (ajj, a) fixed 
for 9m u < 0 and performs the transformation 


( ao \ 


/ (>D \ 

a 

L / rR ; 

— Vo), 

a 

i / /77 ( 


\m/V 2 / 


\m/V 2 / 


(2.17) 


for 9m u >0. It follows that the quantum numbers of the massless state at the singularity 
which were ( 0 , l)i before, now are (M^ 1 ^ J ( 0 ,l)i = (-l,l)o. 

The second case corresponds to the flow from the Nf = 3 theory with three equal bare 
masses mi = m 2 = m 3 = m to the pure gauge theory Nj = 0. Let us discuss in this case the 
quantum numbers of the singularities in more detail. When m = 0, one has two singularities. 
One of them is due to a dyon (1,1) lying in the spinoria.l representation 4 of the flavour group 
Spin(6)=SU(4). A basis for this representation can be taken to be (|0),^|^j|0)) where the 
- 01 , 1 < i < iVy = 3, are the fermionic zero modes carrying flavour indices and one unit of s* 

charge. With this convention, all states with n m = 1 and odd electric charge will lie in the 
same spinoria.l representation 4, while states with n rn = 1 and even n e will lie in the other, 
complex conjugate, spinoria.l representation 4*. The second singularity is due to a n m = 2 
flavour singlet dyon. When n m = 2, the number of fermionic zero modes is doubled and one 


* This would not be the case for k = 1 or k = 2. 

f A senriclassical analysis is valid for studying the flavour symmetry properties of the states becoming 
massless even at strong coupling because these states can be continuously transported to weak coupling 
by varying u. 




Fig. 2 : Shown are three different RG flows from Nj to Nj < Nf , m — > oo, 

Aj\? f —> 0 with A^, A/ = Nf m Nf ~ N f fixed. The upper figure gives the flow from 
f f 

Nf = 2 to Nf = 1, the middle figure the flow Nf = 3 to N'j = 0 and the lower one 
Nf = 2 to N'j = 0. In this latter case which will be studied in detail in Section 

4, we have indicated the quantum numbers of the singularities on both sides of the 
branch cuts. 


can indeed construct SU(4) singlets from the 4 (g) 4* tensor product. The ket associated with 
(1, 2) is then of the form 
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(2.18) 



with 1 < a, b < 2 and a ^ b. Let so be the s = charge of |0). When m > 0, the flavour 
symmetry group is broken from SU(4) to SU(3). We see that the (1,1) singularity will then 
split into one SU(3) singlet of s charge so, (1, l) So , and one SU(3) vector of s charge so + 2 , 
( 1 , l) So + 2 - As for the ( 1 , 2 ) singularity, is must be a SU(3) singlet of s charge so + 3. At the 
expense of shifting the variables ap and a by integer multiples of m/y/2, that is to say choosing 
the cycles ( 71 , 72 ) in (2.5) so that they encircle the poles with non zero residues of the Seiberg- 
Witten differential A an appropriate number of times, we can always choose so = 0 as in Fig. 
2, and the singularities are (1, l)o, (1, 1)2 and (1, 2 ) 3 . This is a natural choice because for 
example the (1, l)o state must remain stable and of finite mass for any u when m —> 00 . This 
choice of s charges will then insure that the solution (ayya) for the Nf = 3 theory will flow 
smoothly towards the solution for the pure gauge theory, without picking any (infinite) shift 
proportional to m. When m = A3/8, we reach a k = 3 superconformal point where the triplet 
(1, 1)2 3 and the singlet ( 1 , 2)3 cross. Using the formula (2.16) with U = + M( 12 ) 3 , 

we see that the singlet ( 1 , 2)3 becomes the dyon (— 1 , 1 )q and the triplet ( 1 , 1)2 becomes the 
quark (1, 0)i. This is exactly what one would expect on physical grounds. Note that we could 
have deduced that the s charge of (1, 2) must be 3 independently of the previous semiclassical 
reasoning by requiring that the elementary quark must have s — 1. 


Finally, the third case we have depicted in Fig. 2 is the flow of the Nf = 2 theory with 
equal bare masses towards the pure gauge theory. This case will be extensively studied in 
Section 4. In particular, we will determine the existence domains of all the BPS states along 
this flow. In the figure, we have also indicated the quantum numbers of the singularities as 
viewed from the upper half u-plane. We have a k = 2 superconformal point when m = A 2 / 2 . 
In the present case, the singlet and doublet collide but do not cross each other and we are 
in the situation of eq. (2.15) with U = M~} = (M( 0 j i) 0 Mn . Note that the upper left 


block of U, acting on (n e ,n m ) is nothing but the matrix —S 



Note also that 


our choice of s charges is not in this case the most natural from the point of view of the RG 
flow. We prefered to make easier the implementation of CP invariance which we will use when 
determining the BPS spectra, by choosing a solution ( 075 , 0 ) satisfying 


ao(u ) = —a]j(u), a{u) = a(u). 


(2.19) 


Then the transformation law of the quantum numbers under CP is simply 

CP 


('M e , Tim)). 


( We,Tim)—s- 


( 2 . 20 ) 
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2 . 3 . The computation of the period integrals 

In the following, we will need explicit expressions for the periods ajj and a, in order to 
compute the curves of marginal stability which determine the existence domains of the stable 
BPS states. A general method, first introduced and used in [8], and also independently in [ 12 ], 
simply consists in expanding the period integrals op and a in terms of the three fundamental 
elliptic integrals, which can be expressed in terms of standard special functions well suited for 
a numerical computation. An efficient way of taking care of the precise definition of the cycles 
(71,72) is to uniformize the cubics ( 2 . 1 ) with the help of the Weirstrafi <p function. For a cubic 
curve in Weierstrafi normal form, i.e. 


V 


2 


3 


e i) j 
i =1 


3 

I= ° 

i =1 


the three fundamental elliptic integrals are 






( 2 . 21 ) 


( 2 . 22 ) 


By convention, we choose the cycle 71 to encircle e3 and e% and the cycle 72 to encircle e\ 
and e2- Of course, we still need to specify how we choose to number the three roots e,, in each 
particular case. Sometimes, if we do not specify the contours, we simply write I\ for f ^ etc. 



Fig. 3 : The definition of the basic cycles 71 and 72 with respect to the roots ei, 
e2 and e3. 


The necessity to introduce the elliptic integral of the third kind I3 is due to the presence 
of poles with non-zero residues in the massive theories. This also explains the failure of the 
standard methods of computation of the periods by means of the Picard-Fuchs equations. For 
the massless theories only I\ and I2 occur which can be reexpressed in terms of hypergeometric 
functions. 
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Let us focus on the case Nt = 2 with mi = m 2 = m. The other cases are discussed in 

jy _2 

Appendix C. The Seiberg-Witten differential A = A m (= m2=m is given by 


A = 


\/2 y dx 


y/2 dx Ay 2 


47t t .2 _ A 2 47t y AH 

X 64 y 2 


X 

V2 dx 
An y 


A ! , A i 

■ x s x ~ a . 


X ~ U + 


A 2 m 2 

4 x 1 M 
I+ 8 J 


(2.23) 


Converting to Weierstrafi normal form of the cubic by y = 2 y, £ = x — ^ we arrive at 


47T 


4 Ai 

-u 1 

3 


1 


21. 


(0 


A 9 z 
—m /. 
2 


2_2 r(*) 


a| 


(2.24) 


7 * 


One sees from (2.23) that A has two poles at (x = —^-,y = Az^mA^) with residues 

[ 2 ]- 

In order to have an explicit formula for ajj and a, we now have to choose the roots ej of 
the polynomial defining the cubic. One constraint comes from the asymptotic behaviour for 
large \u\, which is governed by asymptotic freedom: 


a(u) ~ 



a D (u) ~ 


i 

2n 



(2.25) 


Since a(u) is given by the integral over the cycle 71 surrounding e 2 and e%, see Fig 3, the set 
{e 2 , e%} must be chosen such that the large-u asymptotics of a does not contain a v / 2ulogu/A 2 
term. The remaining root necessarily is e\. But which root in the set {e 2 , 63 } is called 62 , and 
thus is encircled also by 72 , is a matter of convention related to the possibility of shifting aj) 
by an integer multiple of a. A correct choice is the following 


u A| 1 I A§ ~ I A| ~ 

Cl = 6 " 16 + 2V M + T + 2m V M + T _ 2m ’ 



u A| 1 I A| ~ I A 2 ~ 

63 = e " It" 2 V + T + 2m \l U + T~ 2m ■ 


One can then straightforwardly show (see Appendix D) that in the limit m —> 0 the solution 
given by (2.5) and (2.24) converges toward the already known explicit solution of the massless 
Nf = 2 theory [4] which also corresponds to the electric and magnetic quantum numbers at 
the singularities chosen in Fig. 2. It still remains to fix the positions of the cycles relatively 
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to the poles with nonzero residues of A, or, which is equivalent, to fix the s charges at the 
singularities. The solution, consistent with Fig. 2, is given for Am u > 0 by 


alu = 


a D {u) = 


V2 

47T 

V 2 

4:71 


3 U 


T 


2 2 3 


«/| 2) - 2 ir - —m z ir i - 


r( 2) 


A 2 2 r(2) 


A| 


A| 


m 

7! 


(2.27) 


with the precisely given by the formulae 


r(l) 


r(l) 


(ei - ez) l/2 


K(k) 


= 7^—(172 t e i K ( k ) + ( e 3 - e i)E(k)] 


T W ± 

3 (ei - e 3 ) 3 / 2 

4k' 


+ 


1 — c + k' 
1 


K{k) 


1 + k' (1 - c ) 2 - A / 2 


Hi i/(c), 


1 // 

1 + fc' 


(2.28) 


where 


^2 e 2 - e 3 


A; = 


k' 2 = 1 - k 2 = 


c = 


ei - e 3 
c- e 3 
ei - e 3 


z/ c = - 


e2 ~ ei 

e 3 - ei 


1 — c + A/ 
1 — c — k' 


1 — k! 
l + k' 


(2.29) 


and the 1^ obtained form the 1^ by exchanging e\ and e 3 . In (2.28), K, E and III are the 
three standard elliptic integrals of [14] whose integral representations are given in Appendix B. 
For 7 smu < 0 the solution then is obtained from (2.19). In Appendix B we derive eqs. (2.28) 
as well as some useful relations between the elliptic integrals. In Appendix D, we discuss in 
some detail how one can understand the RG flow illustrated in Fig. 2 directly from the explicit 
formulae (2.27). In particular, as m —> oo, m \2 = Ag fixed, we have 


alu) 


*<°)( 


u) 


a D {u) 


i (0) 

'd 


(u) , where ajj[u) = ajj(u) — ea[u) + e 


m 

V2 


(2.30) 


Here and in the following we always let e = sign(Amu). The redefinition ajj —> ajj corresponds 
to rotating the cut originating from the massless quark singularity <r 3 to the right (and changing 
the contour 72 ). This is what one wants since this cut must disappear in the m —> 00 limit in 
order to recover the standard analytic structure of the pure gauge theory. 
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3. The spectra of stable BPS states 


3.1. Decay curves 

In general, a given BPS state does not exist everywhere on the Coulomb branch of the 
moduli space. For each of the massless theories with Nf < 3 [4,3] there is a single curve of 
marginal stability which goes through the singularities and separates the Coulomb branch into 
two regions. In the region outside this curve all semi classically stable states exist, while inside 
the curve only those BPS states exist that are responsible for the singularities, in addition to 
the photon vector multiplet. 

The present cases of hypermultiplets with non-vanishing bare masses are very different. 
Due to the BPS mass formula Mbps = \/2|Z| with 


Z(u) 


nma D {u) 


N f 

— n e a(u ) + 

i— 1 



a BPS state is stable against any decay of the type 

(n e , n m )„^ k X (n' e , n' m ) s ’. + ly. (nf 


(3.1) 


(3.2) 


( k , l G Z) unless this satisfies at the same time the conservation of charges and of the total 
BPS mass: 


n e = kn e + /rig , n rn = kri m + ln m , S{ = ks[ + Is" Z = k Z' + / Z" (3.3) 

and 

\Z\ = \kZ'\ + \lZ"\ (3.4) 

with obvious notations for Z' and Z". If all bare masses mi are equal, clue to the SU(iVj) 
flavour symmetry, only the sum s = ]>T Sj is relevant and needs to be conserved. We see that 
a decay that satisfies the charge conservations (3.3) is possible only if 


§ = C 6 R , (3.5) 

and moreover if it is kinematically possible, i.e. if 

0 < k( < 1 . (3.6) 

For the case of vanishing bare masses, m, = 0 condition (3.5) reduces to 5m = 0 which 
yields a single curve C° on the Coulomb branch independent of the initial state (■ n e ,n m ) Si 
considered. For non-vanishing bare masses however, we have a whole family of possible decay 
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curves. Moreover, a priori, there is a different family of such curves for each BPS state. As an 
example consider a dyon with n m = 1. Then condition (3.5) reads 


n 'm a D - n' e a + Ei s Q | 
a D - n e a + 


0 Sm 


-« - n' m n e )a + 


a D - n e a + JW s 


TYli 

i° l V2 


0 . (3.7) 


For hxed n e and Sj, this is an iVy-parameter family of curves with rational parameters r* = 
(rig ~ ‘n' m n e) /{s'i — n' m Si ). Even though there are some relations between the possible quantum 
numbers n' e and s', n' m (see next subsection) there are still many possible values of r* and 
we expect a multitude of curves of marginal stability on the Coulomb branch of moduli space 
resulting in a rather chaotic situation. Fortunately not all of these curves satisfy the additional 
criterion (3.6). In particular, for the case we will study in detail in the next section, namely 
Nf = 2 with equal bare masses, where one expects a different one-parameter family of curves 
labelled by r = ( n' e — n' m n e )/{s' — n' m s ), s = si + S 2 , for each BPS state, it turned out that 
only one or two such curves in each family are relevant, i.e. satisfy the additional criterion 
(3.6). Hence the set of all relevant curves for all BPS states are nicely described by a single set 
of curves C n E Z, and rather than having a chaotic situation one gets a very clear picture 
of which states exist in which region of the Coulomb branch. It will become clear that this 
organizing scheme should be similarly at work for all other massive theories with Nf < 3. 

One particularly simple case is the decay of states with JJ s i m i — 0 into states with 
JJ s'fxrii = 0. The corresponding decay curves all are given by 9m = 0, i.e. they all 
coincide with the curve C°. We will see that this is quite an important case, and that this 
curve C° still plays a priviledged role, even for non-zero bare masses. 

Note that if we had considered decays into three independent BPS states, (n e , n m ) Si —> 
k x (rig, n' m ) s '. +1 x (n", n'^g" + q x (rig", n(0), s "', we would have two conditions: eq. (3.5) would 
be supplemented by y 6 R, so that such “triple” decays can only occur at the intersection 
points of two curves. Below, when we discuss how to transport a BPS state along a path from 
one region to another, the path can always be chosen so as to avoid such intersection points. 
Hence, triple decays are irrelevant for establishing the existence domains of the BPS states. 
Obviously, “quadruple” and higher decays, if possible at all, are just as irrelevant. 


* This statement is of course meant for the choice of cld and a that gives finite limits under the RG flow as 
m —> oo. For Nf = 2 with equal masses, these are the an and a of eq. (2.30) and JJ SiUii = JJ = 0 
is meant to be s = s' = 0. 
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3.2. Our working hypothesis 

In order to determine the BPS spectra at any point on the Coulomb branch, we will ex¬ 
tensively use the following claim: 

P: At any point of the Coulomb branch of a theory having Nf flavours with bare masses 
rrij, 1 < j < Nf, the set of stable BPS states is included into the set of stable BPS states of 
the m,j = 0 theory at weak coupling. 

Note that the Coulomb branch of the mj = 0 theory is separated into two regions, one con¬ 
taining all the BPS states stable at weak coupling, and the other at strong coupling containing 
a finite subset of the BPS states stable at weak coupling [3,4], One simple consequence of the 
claim (P) is that the set of stable BPS states cannot enlarge when one goes from the Nf to 
the Nf — 1 theory following the RG flow which is what one naturally expects. This is perfectly 
consistent with the spectra determined for zero bare masses in [3,4], Another consequence, 
which plays a prominent role in the present work, is that the possible decay reactions between 
BPS states are then extremely constrained and thus the number of relevant curves of marginal 
stability enormously decreased. This is explained in detail for the Nf = 2 theory with two 
equal bare masses mi = m 2 = m in the next section. 

Let us give a strong argument motivating our fundamental claim (P). It is inspired of ideas 
already discussed in [ 11 ], A crucial fact is that, for any BPS state p = ( n e , n m ) s which does 
not belong to the weak coupling spectrum of the rrij = 0 theory, there always exist some values 
of the rrij for which 

Tfl 

n m aD(uo, rrij) - n e a(uo, rrij ) + = 0 (3.8) 

V2 

at some non-singular uq on the Coulomb branch. This may be proven by noting that for 
sufficiently large rrij |, the point uq , if it exists, lies at large |w| ~ rrij where the formulas for 
ajj and a simplify and thus where (3.8) can be studied very explicitly. The set of curves of 
marginal stability U r£ Q C p (r,mj) a priori relevant for the decays of p = ( n e ,n m ) s , all cross at 
the point uq and form a dense subset of the Coulomb branch. The same properties are true for 
the complementary set U r€ (R_\Q)C p (r, rrij). Now, suppose that p exists in some open set of the 
Coulomb branch. It must then exist at some points on a curve C p (r,mj) with r an irrational 
number. Since the state p cannot decay on such a curve for r irrational, it must also exist at 
the point uq where it is massless, which would contradict the fact that uq is not a singular 
point. Thus p cannot exist as a stable state in any open region for the values of rrij such that 
(3.8) has a solution. Now, by varying the rrij , the curves C p (r,mj) loose their shape and will 
no longer cross at a single point , but we are still insured that p cannot exist on any C p (r,mj ) 
with r irrational. We believe on physical grounds that the fact that p cannot exist in any open 
region of the Coulomb branch means that p simply cannot exist at all as a stable BPS state. 
Finally, note that instead of studying eq. (3.8), we could remark that for arbitrarily small 
rrij |, the spectra of BPS states should be the same as the one for the rrij = 0 theory, which 
proves that an undesirable state like p cannot exist on C p (r, mj) for r irrational and sufficiently 
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small rrij |, and thus for any rrij. However, this reasoning certainly is less rigorous than the 
one based on (3.8). To end this section, let us point out that the claim (P) is also strongly 
supported by the stringy approach used in [ 10 ]. 


4. The case of two flavours with equal masses 


As an illustrative and representative example, in this section we study the BPS spectra 
of the Nj = 2 theory with any mi = m 2 = m real positive bare mass in great detail. The 
analytic structures for small (0 < m < A 2 / 2 ) and large (■ m > A 2 / 2 ) bare masses is displayed 
in the lower Fig. 2, while the explicit solution for the periods ao and a is given by equations 
(2.27) . The three singularities on the Coulomb branch are at points u = crj, o\ < 02 < er 3 , 
such that 


A \ . A| 2 A| 

ai = —— - A 2 m , o - 2 = — 3 - + A 2 m , a 3 = m + — • 

OO O 


(4.1) 


When m = A 2 / 2 , the singularities <72 and <j 3 coincide. For small mass, 0 < m < A 2 / 2 , the 
monodromies M* (see (2.9)) as viewed from the lower half u-plane are given by 


Ml = 

/ 0 1 0\ 

-1 2 0 

, m 2 * = 

/ o 1 0\ 

-12 0 

, M| = 

/ 1 0 0\ 

-2 10 

(4.2) 


\ 1 —1 1 / 


\-l 1 1/ 


\ 0 0 1 / 



while for m > A 2 / 2 , 


/ 0 1 0\ 


/ 2 i 0\ 


n 2 o\ 

-1 2 0 

, M| = 

-10 0 

, m 3 * = 

0 1 0 

V 1 -1 1/ 


V 1 1 1/ 


\0 2 \) 


(4.3) 


The monodromy at the superconformal point, according to (2.11), is given by 



( 

0 

1 

°\ 

M* c = Mf.g = MfM| = 

- 

-1 

0 

0 


V- 

-1 

1 

T —1 


(4.4) 


Before embarking on the detailed derivation of the BPS spectra, let us first discuss what 
one expects from the RG flow arguments. For m — 0 there is a single decay curve, and outside 
this curve all semiclassical states exist, namely all dyons (n, 1) and the W-boson (2, 0) as well 
as the quarks (1,0). The dyons are doublets in one or the other spinor representation of the 
flavour spin(4) group, while the quarks are in the vector representation and the W-boson is 
a singlet. Inside the curve, only the states that can become massless and are responsible for 
the singularities exist, namely the monopole (0,1) and the dyon (—e, 1). As soon as a non¬ 
zero bare mass m is turned on, the s-charge becomes relevant, and certain multiplets split. 
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Semiclassically we still have doublets of dyons (2 n, l)o with even n e and s = 0, while the 
doublets with odd n e split into two singlets (2 n + 1, l)+i and (2 n + 1, l)-i- We have a quark 
doublet (1, 0)i and another quark doublet (l,0)_i, as well as the W-boson (2,0)o- According 
to our claim (P), this is the maximal set of stable BPS states: 

«W= {(2n,l)o 2 , (2n + l) ±1 , (1,0)*?, (2,0) o } • (4.5) 

In the opposite limit, m —> oo, mh .2 = Aq fixed, we expect to flow to the spectrum of the pure 
gauge theory. To describe this limit conveniently, we should move the cut originating from 
the massless quark singularity <73 to the right so that it disappears in the m —> 00 limit. As 
explained at the end of Section 2, we can then choose to work with a and djj instead, which 
are the quantities that flow to and . We call the corresponding quantum numbers n e , 
n m and s : 

n m = n m , n e = n e + en m , s = s + e n m , (4.6) 

where we recall that e = sign (Am u). Hence, in the m —> 00 limit, we expect to have semiclas¬ 
sically the W-boson (h e ,h m )§ = (2, 0)o as well as the dyons (n e ,n m )§ = (2 k, l)o- Converting 
back to the (n e ,n m ) s we always use in this paper, these are ( 2 , 0 )o as well as the dyons 
( 2 n + 1 , l) e : only dyons of odd n e and s = 1 in the upper or s = —1 in the lower half plane 
should survive the RG flow to Nf — 0 in the weak coupling region. All other states must 
disappear. There are two possibilities. A state may simply drop out of the spectrum since 
its BPS mass diverges as m — > 00 as is the case of all states with s ^ 0, i.e s ^ en m . But 
a state can also disappear at a point u (kept fixed ), already at finite m, because it is “hit” 
by its corresponding decay curve which moves outwards as m is increased. We will see below 
that this latter possibility is realised for all dyons (2 n, l)o, n 7 ^ 0 , as well as for the the dyons 
(2 n+ l)_ e , n 7 ^ —1,0. The remaining undesired states, namely (0, l)o, (±1, l)_ e and (l,0)±i 
simply disappear because their BPS masses diverge. 

4.1. Possible decay reactions and decay curves 

In this section, we explain how the claim (P) stated in Section 3.2 drastically restricts the 
number of possible decay reactions and decay curves on which they may occur. We will only 
use the basis of a and ap and the corresponding quantum numbers (n e , n m ) s . Recall that the 
BPS mass is Mgps = \/2| n m aD — n e a + s^=|. As discussed above, the maximal set of BPS 
states then consists of <S max (as well as their antiparticles). These states also constitute the 
semi-classical spectrum. Obviously, if at some point u E M a BPS state decays into two other 
states, the latter must be BPS states and must be contained in this maximal spectrum «S max - 
As in [4] one may also check that matching of the flavour quantum numbers does not give rise 
to any new constraints: the flavour quantum numbers are related to n e and s in such a way 
that they match automatically if n e and s do. The following results do not depend on whether 
the bare mass m is larger or smaller than 4p. For reasons mentioned at the end of Section 3.1, 
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we will only need to consider decays into two types of BPS states: 


(n e ,n m ), -> k x (y.n'roV + l x («"«!!,)«" • 


(4.7) 


The quantum numbers n' m , n'' n always are either 0 or 1, since we can always choose n' m , n" n > 0, 
i.e. a state (n ei —l) s is written as (—1) x (—n e , l)_ s . One has to establish the possible decays 
of each type of BPS state separately. 

As an axample, we present the discussion of all possible decays for the dyons ( 2 n, l)o- The 
corresponding decay curves are 


0 = 5m 


n' m aD - n' e a + s'^= 
an — 2 na 


0 = 5m 


rdQ V2 
an — 2na 


r d = 


n'e - 2 r»C 


(4.8) 


Either n' m = 0 or n' m = 1. Consider first n' m = 0. Then this is the W-boson or a quark. In 
the first case we simply have (2 n, l)o —> k x (2, 0)o + (2n — 2k, l)o- The value of r for this 
decay is oo. In the second case we can have k quarks (1, 0)±i and the dyon (2 n — k, l) s « with 
s" = =F k. Hence \k\ = 1 and we have two possibilities: (2 n, l)o —>► (l,0)±i + (2 n — 1, l)=pi 
and (2n, l)o ^ (—1) x (l,0)±i + (2n + 1, l)±i with r = ±1. Consider now n' m = 1 (and 
n" rl = 1, otherwise we are back to the previous case). Then (2 n, l)o = k x {n' e , 1 ) S ' + / x [v!' e , l) s « 
implies l = 1 — k and ks' + (1 — k)s" = 0. Since |s'|, |s"| = 0 or 1 this implies s' = s" = 0 
(since the cases k — 0 or k — 1 , i.e. I = 0 must be excluded because they do not give 
decays). Hence n' e and rig are even and n = kn' e + (1 — k)n" e . Finally we get (2n, l)o ^ 
k x (2 p, 1) 0 + (1 — k) x (2 q, l)o with n — kp + (1 — k)q with r — oo again. 

We have determined the possible decays for all the other states, namely the W-boson 
(2, 0)o, the dyons (2 n + 1, l)±i and the quarks (1, 0)±i in exactly the same way. The results 
are collected in Table 1. This Table also shows the equations that determine the curves 
corresponding to the decays. 

As already noted above and as will be discussed in detail below, not all curves are relevant. 
It will turn out that the only relevant curves are 


C 


00 


Am 


€Ci Y) 


m 

V2 


= 0 


, ®±7= 

C± : Am-= 0 

ajj — na 


(4.9) 


4.2. Decay curves and BPS spectra for small mass (m < 4p) 

The general picture 

We are now ready to establish the exact existence domains for every BPS state. In this 
subsection, we will discuss the case of small mass, i.e. m < 4p Although it is not too much 
different, it is more convenient to discuss the case of m > 4p separately in the next subsection. 
We will consider each type of BPS state separately. We have seen that for the W-boson and 
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Table 1: 


initial state 


decay products 


W — boson (2, 0)o 
5 m 


r - a °+~ 2 _ Q 


r w = n'Js' 


(I, 0)i + (1, 0)_i r w = 0 

(2p+ 2,1) 0 + (-1) x (2p, 1) 0 r w = oo 

(2 p + 1, l)±i + (-1) x (2 p - 1,1) ±1 r w = ±1 


quarks (1, 0)±i 
_ 


5m 


-“±73 


= 0 


r q = n'm/( s ' =F n 'e) 


(2,0) o + (—1) x (l,0) T i 
(2p, l)o + (—1) x (2p — 1, 1 )=fi 
(2p + 1,l)±i + (—1) x (2p, 1) 0 


r q = 0 

r q = 
r 'q = 


dyons (2n, l)o 


aD—2na 


rd = « - n e n' m )/s' 


k x (2, 0) 0 + (2n - 2k, 1) 0 
(l,0) ±1 + (2 n-l,l) T i 
(-1 ) x (1, 0)±! + (2n + 1,1)±! 
k x (2 p, 1) 0 + (1 - k) x (2g, 1) 0 
with n = kp + (1 — k)q 


r d = cx) 


r d = ±1 
r d = ±1 


r d = cx) 


dyons (2n + 1, l)±i 


k x (2, 0)o T (2n — 2k + 1, l)±i 
(l,0) ±1 + (2n,l) 0 


Ta cl - 

9 ™iR5Tfe = » 2 x(l, 0 ) ±I + ( 2 n-l,l) Tl 


V 2 


rp — 72 e TleTL^j 

d .s' — 


(- 1 ) x (l, 0 ) T i + ( 2 n + 2 ,l ) 0 
(—2) x (1, 0)=|=i + (2n + 3, l)=|=i 


r d = cx 
r d = ±1 
r d = ±1 

r d = Tl 
r d = Tl 


(-1) x (4p - 2 n - 1, l) T i + 2 x (2 p, 1) 0 r d = ±(2n - 2p + 1) 
k x (2p + 1, l)±i + (1 - k) x (2g + 1,1) ±1 
with kp + (1 — /c)g = n r d = oo 


the dyons (2 n, l)o we only need to consider three curves*: r = ±1, oo. For the quarks and 
the dyons (2 n + 1, l)±i we have a priori infinitely many possible decay curves {p G Z), but it 


-k The curve r = 0 for the W-boson or for the quarks is simply 3m a = 0 which is the half-line [ 03 , 00 ). Since 
all states exist in the senriclassical region u —» 00 in the upper and in the lower half plane, every point in 
At can be reached from this region without crossing [ 03 , 00 ). Thus the 3m a = 0 curve is irrelevant. 
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is clear that only a few values of p will correspond to kinematically possible decays: although 
^ G R for all p, only for a few p we will have the additional condition (3.6) that 0 < < 1. 

We have computed all curves numerically using Mathematica. Their exact shape of course 
depends on the value of m/ A2, but is not of much importance. Only their positions relative to 
each other and to the singularities actually matter and there is no qualitative change as long 
as m < y. To check whether a given curve is relevant, i.e. whether 0 < < 1 so that the 

decay is kinematically possible and can really happen, it is enough to proceed as follows: if 
at some point on the curve the decay is kinematically impossible (which can be easily checked 
by computing Z' and Z numerically at this point) one can always transport the BPS state 
through the curve at this point where it cannot decay, and thus the curve is irrelevant. On the 
other hand, to show that a curve is relevant we must make sure that the decay is possible at 
any point of the curve. Since the real-valued function varies smoothly along the curve, it 
is easy to see which real interval is its image and whether it is entirely contained within [0,1]. 

Let us already present the results of the analysis that will be given below. We have 
assembled all the relevant decay curves into Fig. 4 that sketches their relative positions and 
indicates the BPS states that decay across these curves. All curves go through <73, while 
the other intersection point with the real axis depends on the curve: <72, &i and the points 
X2n, n= 1,2,... 

There are several types of states: first, we have the states that become massless at the 
singularities. These are (0, l)o and, clue to the cuts described differently in the two half planes, 
(0, l)o and (—1, l)±i in the upper, and (0, l)o and (1, l)±i in the lower half plane. These states 
exist everywhere (throughout the corresponding half plane). 

Second, we have the other dyons of n e = ±1, the quarks and the W-boson. These states 
decay on curves in the inner, strong coupling region of the Coulomb branch of moduli space: 
The W-boson decays on C°°, the quark (l,0)_i on Cq and the quark (1,0) 1 on the innermost 
curve Cq, while the dyons (e, l)_ e decay on Cq and the dyons (e, l) e on C ( J. 

Third, we have the dyons with \n e \ > 2. As discussed above, among these one must 
distinguish two sorts: those that will survive the RG flow m —» 00 to the pure gauge theory 
and those that do not. The dyons that will survive this RG flow are (2 n + 1, l)i in the upper 
half plane and (2 n + 1, l)_i in the lower half plane. These dyons (n ^ —1, 0) all decay on the 
curve C°° which thus plays a priviledged role. The other dyons, namely (2 n, l)o (n ^ 0) and 
(2 n + 1, l)_i in the upper and (2 n + 1, l)i in the in the lower half plane (n ^ —1, 0) decay on 
curves Cf k , k 7^ 0 (where \2k\ equals \n e \, \n e \ + 1 or \n e \ — 1). There are only two states that 
decay on each of these curves C^ k , k 7^ 0. These curves move more and more outwards as m 
is increased. Also, as \k\ gets bigger (i.e. the \n e \ of the corresponding dyons increase) these 
curves more and more reach out towards the semi classical region. Conversely, as m -> 0, all 
curves flow towards a single curve, say C°°. 


f Here we actually use the fact that the family of curves for a given state is such that the curves do not 
cross each other, except possibly on the real axis. Note that we always consider the parts of a curve in 
the lower and upper half plane separately. 
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Fig. 4 : Shown are a sketch of the relative positions of the relevant decay curves 
for rn < 4^ (for not too large \n e \) as well as the BPS states that decay across these 
curves. Three states do not decay anywhere and still are present in the innermost 
region inside Cq . They are described as (0, l)o and (—1, l)±i in the upper, and as 
(0,1) 0 and (1, l)±i in the lower half plane. Note that, in reality, the angles at which 
the curves meet the real axis at the points are slightly different from what they 
appear to be in the Figure: indeed, the curves Cp fe _ 2 , resp. C£, in the upper half 
plane are the smooth continuations of the curves C+ k , resp. CT 2 , i n the lower half 
plane, in agreement with the monodromy around infinity. 


There are a couple of other points worth mentioning. First remark, that the whole picture 
is compatible with the CP transformation (■ n e ,n m ) s —> (—n e ,n m )_ s under reflection by the 
real w-axis. Second, since all curves go through the singularity 03 , i.e. all existence domains 
touch < 73 , it follows that at this point all BPS states exist. The same is true for the points 
u that lie on the part of the real u line to the right of < 73 . Indeed, as \n e \ is increased, the 
corresponding dyon curves leaving <73 to the right with an ever smaller slope get closer and 
closer to any given point on the real interval (< 73 , 00 ) but never touch it. 

Finally we note that the whole picture is perfectly consistent: if a BPS state decays across 
a given curve, the decay products are also BPS states that must exist in the region considered, 
i.e. on both sides of the curve. Indeed, this is always the case. As an example, consider the 
dyons (2 n, l)o (n > 1). In the upper half plane they decay on the curves C ^ into the dyons 
(2 n — 1, l)i and the quark (1, 0)_i. These dyons (2n — 1, l)i exist everywhere in the upper half 
plane outside C°°, while the quark ( 1 , 0) 1 exists everywhere outside Cq, and in particular in 
the vicinity of the decay curves of (2 n, l)o considered. Moreover, in many cases we perform 
several additional consistency checks to confirm the existence domains we determined. 

Now let us prove that our general picture we just described is indeed true. We begin by 
studying the W-boson. 

The W-boson 



Fig. 5 : Sketch of the possible decay curves of the W-boson (2,0)o for m < 4r- 
The thick curves r w = 1 in the upper half plane and r w = —1 in the lower half 
plane are the only ones that turn out to be relevant: they are C°°. 


From Table 1 we see that the a priori possible decays are into a dyon and an anti-dyon. 
f*w a d H——— 

The three curves 9m —— - — = 0 with r w = ±1, 00 are shown in Fig. 5. 

First consider the curve r w = 00 . On this curve ajj/a is real and varies from —1 to 0 in 
the upper half plane (e > 0) and from 0 to +1 in the lower half plane (e < 0). It is then easy 
to see that the only kinematically allowed decays are (2, 0)o — > (0, l)o + (—1) x (—2, l)o for 
e > 0 and (2, 0)o —> (2, l)o + (—1) x (0, l)o for e < 0. On the r w = +1 curve we similarly 
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find that the kinematically possible decays are (2, 0)o —> (1, l)i + (—1) x (—1, l)i and on the 
r w = —1 curve (2, 0)o —>► (1, l)_i + (—1) x (—1, l)_i for both e > 0 and e < 0. 

Of course, as already mentioned above, such decays can only take place if the final states do 
indeed exist on the curve considered. This will be checked below for the decay on C°°, but for 
the moment we invite the reader to consult the complete Figure 4 to convince herself/himself 
that this is indeed the case. We know that the W-boson exists in region A (connected to the 
semiclassical domain). Also we know that in the m —> 0 limit 02 coincides with o\ and all 
curves become one and the same curve going through o\ = 02 and (73. From [4] we then know 
that the W-boson does not exist inside this single m — 0 curve. By the RG flow to finite m it 
is then clear that it does not exist in the innermost regions D and D'. We will now show that 
the W-boson cannot exist in regions B, B 1 , C or C either. 

Suppose that (2, 0)o exists in region B. Then we can transport it through the cut to region 
B' without crossing any decay curve. There it would be described by (M|. 3 ) _1 (2, 0)o = 
(0, 2)_2- Note that the state (0, 2)_2 is different from (2, 0)o and would have its own decay 
curves. We do not say that (0, 2)_2 would exist in all of region B ', but at least it would exist in 
a region just below the cut that separates B and B'. But such a state with n m = 2 should not 
exist anywhere on M. , hence (2, 0)o cannot exist in B. Suppose now (2, 0)o exists in B'. Then 
transport it through the cut into region B where it is described as M|. 3 (2, 0)o = (0, — 2)_2 which 
again does not exist. Exactly the same argument applies for regions C and C'. Independently 
of the above RG-flow argument, one can see similarly that (2, 0)o cannot exist in D or D' . If it 
would, there would be a (M|) =Fl (2, 0)o = (2, ±4)o in D' or D. We conclude that (2, 0)o cannot 
exist in any of the regions B, B\ C, C' , D and D' . Hence it is the r w = 1 curve for e > 0 and 
the r w = —1 curve for e < 0 that border the existence domain of the W-boson. These are the 
two halves of the C°° curve as defined in (4.9). In this sense, C°° is the only relevant curve for 
this BPS state. 

Note that M2.3 acts on (n e ,n m ) as S G SL(2, Z), i.e. it exchanges n e and n m (up to a 
sign). Thus whenever we have a state with \n e \ > 2 it cannot exist in a region that is bounded 
by part of the cut [01,(72] (as are regions B,B',C and C for the W-boson) since this state 
would be described on the other side of the cut by a \n m \ > 2. This is a very useful fact which 
we will employ much in the following. 

The quark (1, 0)i 

From Table 1 we see that there is an infinity of a priori possible decays on the infinite set 
of curves labelled by r q = — ^, including r q = 00 for p — 0, as indicated in Fig 6. 

One easily checks that for all curves with p ^ 0 the decays are kinematically impossible. 
For p = 0, however, (1, 0)i —> (0, l)o + (—1) x (—1, l)_i is kinematically possible for e > 0 while 
for e < 0 it is (1, 0)i —> (—1) x (0, l)o + (1, l)i. Note that in both cases the quark decays into 


* For typographical convenience we write M*(n e , n m ) s instead of M* 
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Fig. 6: Sketch of the possible decay curves of the quark (1, 0 )i for rn < 4 ^ labelled 
by the values of r q . The thick curve r q = oo is Cq and is the only curve that turns 
out to be relevant. 

the two BPS states that become massless at the singularities a 2 and 03, cf. Fig. 2 , and hence 
indeed exist throughout the corresponding half planes, as we will see below. It is now clear 
that (l,0)i exists everywhere outside the r q = 00 curve which is nothing else than the curve 
Cq as defined in eq. ( 4 . 9 ), while it cannot exist inside, as one sees either from the RG-flow 
argument or otherwise since one would get states (M|) ±:l ( 1, 0)i = (1, q=2)i with \n m \ = 2 that 
do not exist. 

The quark (l, 0 )_i 



Fig. 7: Sketch of the possible decay curves of the quark (l, 0 )_i for m < 4 r 
labelled by the values of r q . The thick curve r q = 00 is C q~ and is the only curve 
that turns out to be relevant. 

Again, there is an infinity of a priori possible decays on the infinite set of curves labelled 
by r q = see Fig. 7 . 

We find again that all decays are kinematically impossible, except for p = 0 on the r q = 00 
curve: (1, 0)_i —> (0, l)o + (-1) x (-1, l)i for e > 0 and (1, 0)_i —> (-1) x (0, l)o + (1, l)-i 
for e < 0 . Note again that the relevant curve goes through <ji and <73 and that the decay is 
precisely into the states that are massless at <ti and <73, cf. Fig 2 . Thus we see that the quark 
( 1 , 0 )_i exists everywhere outside the r q = 00 curve and does not exist inside this curve. The 
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latter fact follows once more from the RG-flow from m = 0 or using (MD ^ 1 which would 
generate a state with \n m \ = 2 that does not exist. Note that here the r q = oo curve is Cq 
which is different from the relevant decay curve C°° for the W-boson (see eq. (4.9)) although 
both curves go through oq and < 73 : Cq lies inside C°°. 

The magnetic monopole (0, l)o 



Fig. 8 : Sketch of the possible decay curves of the magnetic monopole (0, l)o for 
rn < 4r labelled by the values of rq. 

As for all dyons with even n e , there are only three possible decay curves: r = ±1 and 00 . 
For the magnetic monopole, which is massless at <73 these three curves are as shown in Fig. 8 . 

All curves go through <73. Hence any point u can be reached by a path that starts at <73 
and does not cross any curve. Since the monopole certainly exists at (73, we conclude that 
it must exist everywhere. As a consistency check, one can verify that after transporting it 
through any of the cuts it is decribed by quantum numbers that still correspond to allowed 
BPS states that indeed do exist there. 

The dyons (2 n, l)o with n ^ 0 

Let first 2 n > 4. For all these dyons the three curves rq = ±1 and 00 are as sketched in 
Fig. 9. 

In each half plane there are only two curves. The r — 00 curves are the same for all n, 
while the rq = 1 curve in the lower half plane starts at the point X 2 n -2 on the real axis where 
the rq = —1 curve in the upper half plane of the dyon (2 n — 2, l)o had ended. In particular, 
X 2 n increases as n(> 0) is increased. Since the regions B n , B' n , C n and C' n all are bounded by a 
portion of the cut [<7i, < 72 ], by the argument given at the end of the subsection on the W-boson, 
the dyons (2 n, l)o having \n e \ > 2, cannot exist in either of these regions: they can only exist in 
regions A n . The decay on the rq = —1 curve (which is C is (2 n, l)o —> (1, 0)_i + (2 n — 1, l)i 
while the decay on the rq = 1 curve (which is J is (2n, l)o —> (1, 0)i + (2 n — 1, l)_i. 
We will see below that (2 n — 1, l)i, resp. (2 n — 1, l)_i decay on curves that are inside the 
curves C^, resp. C so that the final states of the decays of (2 n, l)o indeed exist. Note 
the following consistency check: if we transport the dyon (2 n, l)o from region A n through the 
interval [x2m x 2 n- 2 ) into the upper half plane, it is described on the other side of the cut as 
M\. 2 . 2 ,{ 2?i, l)o = (2 — 2 n, — l)o = —(2 n — 2, l)o which indeed does exist in the region above 
[X2n, X2n-2\- 
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Fig. 9: Sketch of the possible decay curves of the dyons (2 n, l)o, 2 n > 4, for 
rn < 4r labelled by the values of r f j. The thick curves are those that turn out to 
be relevant. They are C 2n in the upper half plane and C 2n in the lower half plane. 

Now consider the dyon (2, l)o, i.e. n = 1. Then the point X2 n -2 = xq coincides with the 
singularity but this does not change the general argument just given. What is different 
however, is the appearence of an additional curve in the lower half plane, namely r\\ = — 1 
which goes from 02 to <73 and separates region C\ into a region C and an innermost region D. 
The above argument that the dyon cannot exist in region C' n now applies to C. To show that 
(2, l)o cannot exist in D either, it is enough to show that it cannot decay across this r ( \ = —1 
curve separating C and D. The only kinematically possible decay on this curve would be 
(2,l)o ~^ (1, 0 )—1 + (l,l)i- But the quark (l,0)_i does not exist in the neighbourhood on 
either side of this curve as we have shown above. So (2, l)o cannot exist in region D either. 
This also follows from the RG-flow from rn — 0 or by applying when crossing the cut, 
which provides some consistency checks. 

Contrary to the other states considered before, the dyons (2 n, l)o with n > 1 have existence 
domains that are asymmetric with respect to reflection by the real axis, i.e. under complex 
conjugation. However this is not surprising: reflection corresponds to CP and for the states 
considered before, CP only mapped them to their antiparticles which must exist in the same 
domains. On the other hand, CP maps the dyon (2 n, l)o to (—2 n, l)o which is a different BPS 
state. CP then tells us that the dyons (—2 n, l)o exist in domains A n that are the complex 
conjugate of the existence domains A n of the (2 n, l)o- The domains A n are bounded by the 
curves C^ 2n — ^ 2 n the u PP er half plane and by C^ 2n — C 2n in the lower half plane (n > 0). 

Next we will turn to the dyons with odd n e and s = ±1. Since the dyons ( 2 n + 1, l)_i are 


27 



the CP conjugates of the dyons (—2 n — 1, l)+i it is enough to consider the dyons (2 n + 1, l)i 
for all n G Z. For each of these dyons of odd n e , there is an infinity of a priori possible decay 
curves, one for every odd r and one for r — oo. 

The dyon (—1, l)i 


i 



Fig. 10: Sketch of the possible decay curves of the dyon (—1,1) i for m < ^re¬ 
labelled by the values of r ( \. The thick curve (Cj~) is the one that turns out to be 
relevant. 


All curves go through < 73 , and the curves in the upper half plane all end at cr\ where this 
dyon is massless. The curves are shown in Fig. 10. 

Since the dyon (— 1 , l)i exists at a\ when approached from the upper half plane, we can 
always transport it to any point in this half plane without crossing any curve, and thus we know 
that it exists everywhere in the upper half plane. Let us now focus on the lower half plane. The 
curves with r ( \ = — 2p — 1 correspond to the decays (—1, l)i —> 2x (2 p, l)o + (—1) x (4p+l, l)_p 
Clearly, if \p\ is too large, the final states will be too massive and the decay is impossible. We 
find that on the r\] = —1 curve (p = 0 ) the decay (— 1 , l)i —> 2x ( 0 , l)o+(— 1 ) x ( 1 , l)_i is indeed 
possible, i.e. (— 1 , l)i decays into the states that are massless at (73 and <71 in the lower half 
plane. We also find that no decays are kinematically possible on the r ( \ = —3, —5, ... curves. In 
any case, (— 1 , l)i cannot exist in any of the regions in the lower half plane between two curves 
rd = —2 p — 1 and rd = —2 p — 3 (p = 0, 1, 2, ...) because these regions touch part of the cut 
[(7 1 , < 72 ] and hence the existence of (—1, l)i would imply the existence of M|. 3 (—1, l)i = ( 1 , 1)3 
just above the cut. However, such a state has s = 3 and cannot exist. Similarly, (—1, l)i 
cannot exist between any of the curves rd = 2p + 1 and rd = 2p + 3 (p — 0, 1, 2, ...) in the 
lower half plane because they touch part of the cut [( 72 ,( 73 ], and we would similarly conclude 
that a state M 3 (—1,1) 1 = (—1, 3) 1 would exist just above the cut. We conclude that the dyon 
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(— 1 , l)i exists everywhere except in the region of the lower half plane that is bounded by the 
cut [< 7 i, 0 - 3 ] and the r ( \ — —1 curve which is nothing else than Cq (or actually its part in the 
lower half plane). 

As discussed above, it follows from CP that the dyon ( 1 , l)_i exists everywhere outside 
the mirror image of the region just described, which is bounded by [<Ji, 03 ] and the part in the 
upper half plane of the curve Cq. 

The dyon (1, l)i 



Fig. 11: Sketch of the possible decay curves of the dyon (1, l)i for m < 4r labelled 
by the values of r ( \ . The thick curve (Cq) is the one that turns out to be relevant. 

All curves go through < 73 , and the curves in the lower half plane all start at 02 where this 
dyon is massless. They are shown in Fig. 11. 

Again, any point in the lower half plane can be reached from 02 without crossing any curve, 
and hence (1, l)i exists everywhere in the lower half plane. What about the upper half plane? 
It turns out that there the only curve on which a decay is kinematically possible is the r c j = 1 
curve with the decay ( 1 , l)i —> 2 x ( 0 , l)o + (— 1 ) x (— 1 , l)_i with the final states being again 
those BPS states that are massless at 03 and 02 in the upper half plane. It is then clear that 
this decay must happen as (1, l)i cannot exist inside the r ( \ = 1 curve. This can be seen either 
from the RG-flow argument from m = 0, or else since the existence of (1,1) 1 in this region 
would imply the existence of (M|) _1 ( 1 , l)i = ( 1 ,3)i below the cut [ 02 , 03 ] which is excluded. 
We conclude that the dyon (1,1) 1 exists everywhere outside the region bounded by [ 02 , 03 ] 
and the r ( \ = 1 curve in the upper half plane which is nothing else than the corresponding part 
of Cq . By CP we see that the dyon (—1, l)_i exists everywhere outside a region bounded by 
[ 02 , 03 ] and the part of Cq in the lower half plane. 

The dyons (2 n + 1, l)i with n > 0 
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Fig. 12: Sketch of the possible decay curves of the dyons (2 n + 1, l)i for m < ^re¬ 
labelled by the values of r\\. The thick curves are those that turn out to be relevant. 

They are C°° in the upper half plane and in the lower half plane. 

Among the infinity of a priori possible decay curves, we find that the decays are kinemati¬ 
cally possible only on the following ones: in the upper half plane on r ( \ = oo, r ( \ = 2n + 1 and 
r,\ = 2 n — 1 and in the lower half plane on rq = oo, r\\ = 2n — 1 and r ( \ = 1. These curves are 
shown in Fig. 12. 

The curves shown in Fig. 12 are generic, except for 2n + 1 = 5 where the r\] = 2n — 1 = 3 
curve in the upper half plane ends at oq, and for 2n + 1 = 3 where there is no r ( \ = 2n — 1 = 1 
curve in the upper half plane and where the r ( \ = 2n — 1 curve in the lower half plane coincides 
with the r ( ] = 1 curve and starts at xq = <7i. 

There are two kinematically possible decays on the r ( \ = 1 curve: (2 n + 1, l)i —> (1, 0)i + 
(2 n, l)o and (2 n + l,l)i —>► 2 x (1,0)i + (2 n — 1, l)_i. But this curve is and hence also 
the decay curve of (2 n, l)o (where the latter decays into (1,0)i + (2 n — 1, l)-i), so that the 
first of the two decays cannot take place (or actually is identical to the second). On the other 
hand, we will see below that (2 n — 1, l)_i exists in the lower half plane everywhere outside 
the curve C°° which is well inside (n > 0), so that the second of the two decays can indeed 
take place. In the upper half plane, on the r ( \ = oo curve (which is C°° ) we can have the decay 
(2 n + 1, l)i —> (n + 1) x (1, l)i + (— n) x (—1, l)i into final states that do exist there. The 
existence of (2 n + 1, l)i in regions B, C , D, B\ C and D' is ruled out by the by now familiar 
arguments of transporting the state through the cuts [ 01 , < 72 ] or [< 72 ,( 73 ] and thus generating 
states (M|. 3 ) ±:L (2n + 1, l)i or (M|) ±:L (2n + 1, l)i that would have |n m | > 2. 

The dyons (2 n + 1, l)i with n < —1 

The decay curves on which the decays are kinematically possible look qualitatively very 
much like those of Fig. 12. The only curve that starts in the lower half plane to the left of aq 
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is the curve C 2n , 2 (which is now rvi = —1 and replaces the r\\ = 1 curve of Fig. 12), while in 
the upper half plane we still have the same r — oo curve C °°. By the same arguments, using 
(Mf.s)* 1 and (M|) ±:L , we can show that the dyons under consideration only exist outside these 
two curves. The decay on C°° is e.g. (2 n + 1, l)i —► (n + 1) x (1, l)i + (— n) x (—1, l)i as 
before, while on C 2n+2 it is (2 n + 1, l)i —»• 2 x (—1, 0)i + (2 n + 3, l)_i which do exist there. 

The dyons (2n + 1, l)_i 

These dyons are the CP conjugate states of (—2 n — 1, l)i. Hence, in the lower half plane 
they always exist everywhere outside the curve C°°, while in the upper half plane they only 
exist outside curves C 2n if n > 0 and only outside curves C 2n+2 if n < — 1. As a consistency 
check, consider transporting a dyon (2 n+ 1, l)i with n > 0 through the cut (—oo, 01 ] (see Fig. 
12) into the region below the cut. Then this state is described there by (M^. 2 . 3 )~ 1 (2ni, l)i = 
— (2 n + 3, l)_i. This state indeed exists in this region because it is outside C°° (in the lower 
half plane). 

4.3. Decay curves and BPS spectra for large mass (m > 4r) 

The general picture 

As m is increased from m < y to m > y one goes through rri — 4 ^ where the singularities 
02 and 03 coincide and where new states become massless. At m > 4 ^, the quantum numbers 
of the massless states at these two singularities have changed. One must however keep in mind 
that this is a somewhat “local” effect, where locality here refers to the distance on the Coulomb 
branch: these rearrangements do not much affect a(u) or ao(u) for u far away, i.e. such that 
\u — 02 1 3 > 1 02 — 03 1 . What does this mean for the relevant decay curves and existence domains 
of the BPS states of Fig. 4 ? We expect, and verify below, that all what can happen is the 
following: as m = 4 jr arid 02 = 03, the curve Cq has shrunk to a point. Then as m is increased 
beyond 4 jr and 03 moves off to the right, all curves (except Cq") remain “attached” to the real 
axis “at the left” at points oq, X2 n (that, of course, do move as m is varied), while at their 
right some curves remain attached at 02 while others are attached at 03. It turns out that all 
curves C 2n , n G Z are attached to 02, and all curves C 2n with n < 0 for e > 0 and n > 0 for 
e < 0 are attached to 03. But there appears also a new family of relevant curves, attached to 
02 and 03: these are the curves C 2n with n > 0 for e > 0 and n < 0 for e < 0. All this is shown 
in Fig. 13 . The only states existing inside the Cq" curve are the three dyons (e, 1 ) € , (—e, l) e , 
(—e, l)_ e as well as the quark ( 1 , 0 ) 1 and the monopole ( 0 , l)o- The dyon (—e, l)_ e and the 
monopole ( 0 , l)o however decay on the curve Cq between 02 and 03. So, as always, the only 
BPS states that exist everywhere are the three states that are responsible for the singularities, 
namely the two dyons (±e, l) e and the quark (1, 0)i. 

Note that, as m —■> 00 , the special dyons (—e, l) e and (e, l) e have finite masses as do all 
dyons (2 n + 1, l) e as well as the W-boson. These states thus survive the RG flow to the pure 
gauge theory. These surviving states all decay on one and the same curve, C°° which one may 
consider as fixed under the flow (except for (—e, l) e and (e, l) e , of course). After transforming 
to the natural quantum numbers ( n e ,n m )§ of the pure gauge theory, see eq. (4.6), we see that 
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Fig. 13: Shown are a sketch of the relative positions of the relevant decay curves 
for rri > 4^ as well as the BPS states that decay across these curves (for not too 
large \n e \). The BPS states that decay across C°° in the upper half plane are the W- 
boson (2, 0)o as well as all dyons (2 n + 1 , l)i, n ^ — 1 , 0, while the states decaying 
across C°° in the lower half plane are again the W-boson (2, 0)o as well as all dyons 
(2 n + 1, l)_i, n ^ —1,0. The only states existing inside the Cq curve are the 
three dyons (e, l) e , (—e, l) e , (—e, l)_ e as well as the quark (1, 0)i and the monopole 
(0, l)o- 
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these states all have s = 0 and precisely constitute the spectrum of the pure gauge theory as 
established in [3]. On the other hand, the states that do not survive this flow either simply 
disappear from the spectrum because their masses diverge, as is the case of the quarks or the 
special dyons (±1, l)_ e or (0, l)o inside C°°, or because they are “hit” by their decay curve 
that moves outwards as m —» oo. Indeed, any point u which for some m still is in the existence 
domain of a given BPS state will be hit by the corresponding curve for large enough m. Note 
that again, as for rn < 4p, all BPS states exist at <73 and in a narraow wedge extending to 
the right of it, in particular on the real interval [< 73 , 00 ). Of course, this does not prevent any 
real point u to be hit by the curves: As m —> 00 we also have (73 —> 00 and any point u will 
eventually end up inside the curves. 

Now let us briefly comment on the discrepancy with the published version of [9]. The 
authors of [9], for Nf = 2 with equal bare masses, consider a real point u to the right of 
the singularity (73 and write that a certain dyon does no longer exist there for a certain mass 
m > 4r- In the light of what we have said so far, it is clear that this cannot be true. This fact 
is even more obvious since, if it were true, we would need to have decay curves crossing the 
real u axis to the right of <73 where there are no cuts - a situation that cannot occur as one 
can easily see, also without our detailed analysis. The authors of [9] have checked this point 
again, and told us that they actually agree with our result, the discrepancy being only due to 
some error when writing up their paper. 

In the remainder of this subsection, we will again discuss the decay curves for each type 
of BPS state, thus proving Fig. 13 to be correct. 

The dyons (2 n, l)o 

Recall that decays can only happen on the curves (4.8) for r = ±1 and 00 . But numerically 
we find that the ryj = 00 curve has disappeared. Actually, as m approaches from below, 
one can see how this curve (which goes through < 73 , see Figs. 8 to 9) shrinks to a point, and 
its absence at m > is perfectly compatible with a smooth RG flow as m is increased. 

For a generic dyon (2 n, l)o with n > 2 one finds for the r f j = ±1 curves the situation 
depicted in Fig. 14 on the left. The r ( \ = 1 curve (for e > 0 and e < 0) is the curve C ^ 
while the r ( \ = —1 curve (e > 0) is C^. All curves turn out to be relevant with the following 
kinematically possible decays: (2n, l)o —> (l,0)_i + (2 n — l,l)i on r& = —1, (2n, l)o —>► 
(—1) x (1, 0)i + (2n+ 1, l)i on r\\ = 1 for e > 0 and (2 n, l)o —> (1, 0)i + (2 n — 1, l)_i on rq = 1 
for e < 0. To show that (2 n, l)o cannot exist in regions C n , B n and B' n we use the same type 
of argument as before. The monodromy M2.3 around <72 and <73 is the same as for small mass, 
and just as before, transporting the dyon (2 n, l)o from C n to B' n through the cut [(71,(72] or 
vice versa would result in states with |n m | = 2 n that cannot exist. To prove the non-existence 
of (2 n, l)o in B n , transport it to B' n where it is described as (M|) _1 (2n, l)o = (2 n — 2,1)_2 
which has s = —2 and cannot exist. We conclude that (2 n, l)o only exists in regions A n 
outside these curves. The only difference for n = 1, i.e. for the dyon (2, 1 )q is that the point 
%2n-2 = xq coincides with oq, but this does not change the preceeding discussion. 

The dyons (2 n, l)o with n < 0 are the CP conjugates of the dyons (—2 n, l)o with —2 n > 0, 
and thus exist in the complex conjugate domains A n = A_ n . 
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Fig. 14: In the left figure we show a sketch of the relative positions of the decay 
curves of the dyons (2 n, l)o (n > 1) for m > 4p. All curves shown are relevant. 

The hgure on the right shows a sketch of the relative positions of the decay curves 
of the monopole (0, l)o- The thick curves are the relevant ones. 

The situation is slightly different for the monopole (0, l)o, as shown in Fig. 14 on the 
right. No decays are possible on the r = —1 curve while on the r —l curve one has (0, l)o —>► 
(— 1) x (1, 0) i + (1, l)i for e > 0 and (0, l)o ~^ (1, 0) i + (—1, l)-i for e < 0. The monopole 
cannot exist in B or B' since this would lead to (M|) ±1 (0, l)o = (±2, 1)±2 on the other side 
of [o 2 ,C 3 ] which, having |s| = 2, is excluded. 

The dyons (2 n + 1, l)i 




Fig. 15: Shown is a sketch of the relative positions, for m > 4r, of the decay 
curves of the dyon (1, l)i (left) and of the dyon (—1, l)i (right). The thick curves 
are the relevant ones. 

First, look at the dyon (1, l)i for which the curves are as shown in the left Fig. 15. All 
curves go through <72 where this dyon is massless when approached from e > 0, and we see 
that it must exist in the whole upper half plane. It cannot exist in region B' since otherwise 
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this would imply the existence of M|( 1, l)i = (3, 1)3 above the cut. Note that conversely, the 
existence of (1, l)i in B only implies the existence of (M|) _1 ( 1, l)i = (—1, l)_i below the cut 
which is perfectly consistent (see Fig 13). The decay across the r d = 1 curve for e < 0 (Cq ) is 
( 1 , l)i — > 2 x (l, 0 )i + (— 1 , l)-i into final states that are massless in the lower half plane at 
<73 and ( 72 . 

Next, consider the dyon (—1, l)i with curves shown in the right Fig. 15. All curves in the 
upper half plane go through cri where this dyon is massless. Hence it exists everywhere in the 
upper half plane. It cannot exist in B or in C since M|. 3 (—1, l)i = (1, 1)3 and M|(—1, l)i = 
(1, 1 ) 3 . The decays in the lower half plane are (—1, l)i —> (—1) x (1, 0)_i + (0, l)o on r d = —1 
(Cq') and (- 1 , l)i-> 2 x ( 1 , 0 )i + (-3, l)_i on r d = 1 (CZ 2 )- 




Fig. 16: Shown is a sketch, for m > 4jr, of the relative positions of the decay 
curves of a generic dyon (2 n + 1 , l)i with n > 1 in the left figure and n < —2 in 
the right figure. We have not shown the curves inside the r — 00 curves. The thick 
curves are the relevant ones. 

Now, for a generic dyon (2n+1, l)i with n > 1 one has the situation of the left Fig. 16. Of 
course, for n = 1, i.e. for the dyon (3,1) 1, the r — 1 curve for e < 0 starts at xo = a 1, but this 
makes not much difference. Decays are only possible on the r — 00 curve and the r — 1 curve 
for e < 0. By the same arguments as above it is clear (for n > 1) that (2 n + 1, l)i cannot exist 
in regions B n or C n . The decays are (2 n + 1, l)i —> (n + 1) x (1, l)i + (— n) x (—1, l)i on the 
r — 00 curve which is C°° in the upper half plane, and (2 n + 1 , l)i —> 2 x ( 1 , 0)i + ( 2 n — 1 , l)_i 
on the r — 1 curve which is C 2n in the lower half plane. 

Next, for the dyons (2 n + 1, l)i with n < —2, we have the situation of the right Fig 16. 
The kinematically possible decays are: (2 n + 1, l)i ->2x (1, 0)i + (2 n — 1, l)_i on the r d = 1 
curve for e < 0 (C 2n )^ (2 n + 1, l)i —> (—2) x (1, 0)_i + (2 n + 3, l)_i on the r d = —1 curve for 
e < 0 (C^ +2 ), and ( 2 n + 1 , l)i —> (—n) x (— 1 , l)i + (n + 1 ) x ( 1 , l)i on the r d = 00 curve 
for e > 0 (C 00 ). By the (M|. 3 ) ±:L and M| arguments one sees that (2 n + 1, l)i with n < —2 
cannot exist in regions B n: C n and D n . 

The dyons (2n+ 1, l)_i exist in the complex conjugate domains of where the (—2 n — 1, l)i 
exist. All this is shown in Fig. 13. 
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The quarks (1, 0)±i 




Fig. 17: Shown is a sketch, for m > 4jr, of the relative positions of the decay 
curves of the quark (1,0)i (left) and of the quark (l,0)_i (right). For the latter, 
we have not shown the curves r > 0 inside the r — oo curve. The thick curves are 
the relevant ones. 


First, the quark (1, 0)i is massless at < 73 . All curves go through < 73 , see Fig. 17, and we 
conclude that this quark exists everywhere! 

For the other quark (l,0)_i, the curves are shown in the right Fig. 17. Decays are only 
possible on the r — 00 curve (Cq") where they are ( 1 , 0)_i —> (0, l)o + (— 1 ) x (— 1 , l)i for e > 0 
and (l,0)_i —> (—1) x (0, l)o + (1, l)_i for e < 0. The quark (l,0)_i cannot exist inside the 
r — 00 curve since again this would imply the existence of a state (M|. 3 ) ±1 (l, 0)_i = (0, =Fl )_2 
that does not exist. 

The W-boson (2, 0)o 

For the W-boson, there is no r — 00 curve, and the r = ±1 curves exist only for e = ±1, 
where they actually are C°°. By the (M|. 3 ) ±:L argument, the W-boson cannot exist inside C°°. 
The decays are (2, 0)o —» (1, l) e + (—1) x (—1, l) e . 
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5. Physical discussion : the N = 2 superconformal fixed points 

In the light of our preceding results, we now discuss the physics of the superconformal 
points. 

At a point -u* on the Coulomb branch where mutually non-local dyons become massless, 
the low energy theory is believed to be an interacting superconformal held theory [ 6 ]. In 
the theories we studied in this paper, such a point occurs when two singularities coincide 
[7]. Near w* which always lies in a region where the low energy theory is strongly coupled 
(independently of the choice of variables), the masses of the particles becoming massless at w* 
are much smaller than the masses of all the other particles (which are of order A) and set a 
new mass scale M. This implies that one can give an intrinsic meaning to the superconformal 
held theory (SOFT) independently of its embedding in the original non-abelian gauge theory, 
by letting A —> oo. Actually, the same N = 2 SOFT can be embedded in different N = 2 
non-abelian gauge theories. Argyres and Douglas were able to show, in a particular but generic 
case that corresponds in our language of Section 2.1 to a k — 1 SOFT, that the low energy 
coupling does not depend on the separation of scales A/M. This is a strong indication that the 
coupling does not run between the scales M and A, and thus very convincing evidence that the 
theory indeed has conformal invariance. One of the most striking properties of these SCFTs 
is that they do not contain any gauge boson that could contribute to the (3 function with a 
minus sign. We have seen this in great detail in Section 4 for the k = 2 superconformal point 
appearing in the Nf = 2 theory, since the only spin one particle at this point has quantum 
numbers (2, 0)o and a mass A 2/\/2 ^ 0. It is easy to realize that this is also the case at any 
superconformal point appearing in the 1 < Nf < 3 theories. To explain that the [3 function 
can nevertheless be zero, Argyres and Douglas [ 6 ] suggested a simple ansatz which consists 
in computing the perturbative contribution of each hypermultiplet which occurs in the SOFT 
separately. This is done by using duality to go to a formulation of the theory where the given 
hypermultiplet is described locally, and then to apply the inverse duality transformation to 
the contribution to the beta function. The individual terms are then simply added to obtain 
the total (3 function: 

E("-» t -^) 2 - ad 

where each state (n e ,n m ) s in the sum should be counted with its correct multiplicity d. Of 
course, it is by no means obvious why this should give the correct answer. In this equation, 
r = 9/7 T + 8in/g 2 is the (generalized) coupling of a theory containing the hypermultiplets 
(n e , n m ) over which the sum is done. This is not the low energy coupling of the original gauge 
theory which, to avoid confusion, we will hereafter denote r e fj = dap/da . 

To discuss the validity of this ansatz (5.1) remark that it has three immediate and impor¬ 
tant consequences noted in [ 6 ]. The first one is that the 6 angle does run, due to the contribu¬ 
tion of magnetically charged states. This is quite surprising in a U(l) theory, and this effect 


* An early and somewhat related discussion can be found in ref. 22. 
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may well remain qualitatively valid and have interesting consequences in non-supersymmetric 
theories as well, perhaps by shading a new light on the strong CP problem. However, we will 
have nothing new to say about this in the following. The second consequence of (5.1) is that 
we have an RG fixed point with a fixed point coupling r* satisfying 



r* — n 


2 


0 . 


(5.2) 


The third consequence is that this fixed point is IR stable, the slope u> of the /3 function being 
positive, with (5.1) giving 


1 <9(t — T*) (X ^ 

u = lie hm-—-=- hmr* > 0. (5.3) 

T—>75 * r — T* O log /I 7 T 

The positivity of u ; is also required by unitarity [15]. In the following, we will argue that (5.2) 
is correct, but that (5.3) is wrong. This was already suspected in [6], since (5.3) gives irrational 
values for c o. We will see that the correct value for u is indeed a rational number. 

Before presenting our discussion, it is necessary to stress the following: although the 
effective coupling r e fj(u) and the running coupling r u (/i) of the microscopic theory at fixed u 
and scale /i are different physical quantities, they are nevertheless related. The easiest way to 
understand this is to first choose u near a singular point uq where only locally related states 
become massless, a familiar case. At fixed u, the particle spectrum of the theory consists 
of the BPS states becoming massless at uq and having masses M ~ t(u — uq ) where t is a 
constant, and of other states of mass ~ A or higher. The dependence of M on u simply means 
that u is a good local coordinate near uq, or, in terms more appropriate for the discussion to 
come, that u has dimension one. This is equivalent to the fact that the low energy theory 
is free massless super-QED. The effective coupling r e ff for a given u then corresponds to the 
microscopic running coupling t u (ii) at a scale n ~ M or lower. Between the scales M and A, 
r u (/r) is given reliably, after a duality transformation which renders the theory weakly coupled 
in the IR, by a one loop calculation in ordinary N = 2 super-QED coupled only with the states 
of masses ~ M. One then identifies r e ff(u) with r w (yU = M) to End the asymptotic behaviour 
of r e ff near the singularity [1,2]. Similarly, using the same argument backwards, and repeating 
this reasoning near a superconformal point u* where r e fj is known via our explicit formulas 
(e.g. (2.27) for the k = 2 SOFT), one can deduce the form of r M (/x). Hence the effective 
coupling at w* must coincide with the fixed point coupling r*: 

TeffCu*) = T* . (5.4) 

Furthermore, we will see below that r e fj has the following asymptotic behaviour, 

'Tefr(w) = T* + C(u — u*y + o((w — w*) 7 ). (5.5) 

Moreover, it follows from eq. (2.13) that at u the masses of the particles becoming massless 
at w* are M ~ t{u — w*) 1 /®, and that a is the anomalous dimension of u. We deduce that 
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r — r* ~ p° 7 and thus that the slope u> of the /3-function is given by 


u = a'f. (5.6) 

Let us first exploit the fact that r e ff(u*) = t* to discuss (5.2). Naively, at a superconformal 
point where two singularities coincide, one has precisely two massless states. For instance, at 
the k — 3 superconformal point of the Nj = 3 theory depicted in Fig. 2, one would expect to 
have the states (1, 1)2 and (1, 2 ) 3 . However, with these states, (5.2) does not give the correct 
answer for the fixed point coupling, which in this case is r* = e * 71 "/ 3 = r e ff(u*). The solution 
to this puzzle is not to give up (5.2), but to realize that there are other states which become 
massless at the superconformal points. One way to find them is to look for massless states in 
the maximal set of BPS states of the Nf = 3 theory. This is correct because one can show, 
along the lines of Section 4, that all the states belonging to the maximal set must exist at 
the superconformal point. Thus we deduce that we have six massless states, three triplets 

( 1 . 1 ) 2 , (0,l)i and (l,0)i, and three singlets ( 1 , 2 ) 3 , (—1, l)o and ( 2 , 1 ) 3 . With these states 

(5.2) gives the right answer. One interesting point to note is the following: the states being 
massless at the potentially coinciding singularities for small masses (m < A 3 / 8 ), i.e. (1, 1)2 and 
( 1 , 2 ) 3 , and the states becoming massless for large masses (m > A 3 / 8 ), i.e. (— 1 , l)o and ( 1 , 0) 1 , 
are four distinct states and not simply analytic continuations of each other. The mechanism 
which allows to change the quantum numbers at the singularities in the asymptotically free 
theories, a phenomena which was known to occur since the work of Seiberg and Witten [1,2], 
is thus fully elucidated: at a point where two singularities collide, new states become massless 
“accidentally” and, when the singularities split again, these new states can take over the roles 
of the original massless states which were responsible for the singularities before the collision. 

Let us illustrate this mechanism in more detail in the case of the Nf = 2 theory studied 
in Section 4. Figure 18 represents an enlargement of Figs. 4 and 13 near the singularities. 
We explicitly indicate the existing states in the different regions surrounded by the curves of 
marginal stability. Since a D = a — = 0 at the superconformal point, it is clear from Figure 

18 that the states becoming massless at this point are the two singlets (—1, l)_i and (1, l)i, 
and the two doublets (0, l)o and (1,0) i- Equation (5.2) then gives the correct fixed point 
coupling, r* = i. When m < A 2 / 2 , the states responsible for the singularities are (—1,1 )—1 
and (0,l)o, and exist everywhere on the Coulomb branch. The states (1,1) 1 and (1,0) 1 are 
also very special, since they exist everywhere except inside the innermost curve Cq . This is 
a nice test of the spectra derived in Section 4, since (1, l)i and (1, 0) 1 , being massless at the 
superconformal point, must exist everywhere except possibly inside curves of marginal stability 
shrinking to zero when such a point is approached. This is indeed what happens for the curve 
Cq. When m > A 2 / 2 , the roles of the states (—1, l)_i, (0, l)o and (1, l)i, (1, 0)i are exchanged 
but the picture still is perfectly coherent because (—1, l)_i and (0, l)o again exist everywhere 
outside Cq. 

Let us mention that at the k — 1 SCFT appearing in the Nf — 1 theory, we have three 
massless states which correspond exactly to the three states which were identified in [ 6 ] at the 
Z 3 point of the SU(3) moduli space. This was expected since the two SCFTs are believed [7] 
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(- 1 , 1 ) 


( 1 , 1 ), 


( 1 , 0 ), 




Fig. 18: The curves of marginal stability and stable BPS states in the upper half 
u -plane near the singularities in the Nf — 2 theory with mi = m 2 = m. The upper 
configuration corresponds to m < A 2/2 and the lower configuration to m > A 2 / 2 . 


to be the same. It would be interesting to check whether such a correspondence between states 
still holds for the k = 2 and k — 3 SOFT when these theories are embedded in the SU(4) or 
SO (8) pure gauge theories [16]. 
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There is still one unsatisfying point in the above discussion. Though we saw that (5.2) 
gives the correct fixed point coupling, its application requires the knowledge of the set of all 
the massless states, which from the above discussion seems to depend crucially on the theory in 
which the SCFT is embedded. This is not natural since r* is, at least modulo a duality trans¬ 
formation, a characteristic of the SCFT. Actually, the massless spectra can be deduced from 
the local data only. Indeed, as the original massless states exist everywhere on the Coulomb 
branch and cannot decay, one can produce new massless states by encircling the superconfor- 
mal singular point. The massless spectra must then be invariant under the monodromy M sc 
at the superconformal point. This invariance only generates a finite number of new states 
because M sc is of finite order. Using this invariance, one can deduce straightforwardly the 
massless BPS spectra uniquely from the quantum numbers of the colliding singularities. For 
instance, by applying (4.4) on the states (—1, l)_i and (0, l)o, one generates exactly the two 
missing states (l,l)i and ( 1 , 0 )i. 

It is very tempting to conjecture that the matrix M sc always corresponds to an exact 
quantum duality symmetry of the SCFT. In addition to the invariance of the spectrum, this is 
supported by the (related) fact that the fixed point couplings satisfy M sc ■ r* = r* and by the 
fact that the symmetry M sc can be realized as a global symmetry e.g. as the Z 3 symmetry of 
the k — 1 SCFT when embedded in the SU(3) pure gauge theory. Note that the surprising idea 
that global symmetries can be related to non-trivial duality transformations was first realized 
and used in [3] in the case of the Z 2 symmetry acting on the moduli space of the SU(2) pure 
gauge theory. We find a very explicit realization of this fact with the superconformal points 
studied in this paper. Very recently, this idea has also been used very nicely to study S duality 
in some finite gauge theories [17]. We strongly feel that this is a very interesting point of view, 
since it relates a priori purely quantum symmetries to more conventional symmetries which 
have a classical origin. 

Let us now turn to the computation of c 0 . We will also obtain a, recovering the results of 
[7]. To compute these two critical exponents, it is enough to find the asymptotic expansions 
of a 75 and a near -u* up to the second order. An easy way to deduce the result is the following: 
diagonalizing the monodromy matrix M sc , its eigenvalues are x± = exp(±^(l + k )), which 
shows that the leading power in the expansion, corresponding to 1/ck, see (2.13), must be of 
the form 

1 . 1 + k 

- = ±—— + n, ne Z. 

a 12 


Imposing a > 1 which is a necessary condition for a unitary SCFT [7], and a < 2 which 
amounts to saying that the operator u is relevant, the only possibility is that a is related to 
the eigenvalue X- with n — 1 , 


a = 


12 

11 - k 


(5.7) 


The other eigenvalue x+ is related to the subleading power 1/(3 > 1/a in the expansion by 


1 

1 



13 +A; 
12 
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(5.8) 



The formulae (5.7) and (5.8) can be checked using the explicit solutions given in (2.27) for 
Nf = 2 or those in Appendix C for the other cases. By using 


CLd(u ) =CLd(u*) + C£)(u — W*) 1 ^ + c D (u — M*) 1///3 + o((u — M*) 1///3 ) 
a(w) =a(«*) + c(u — -u*) 1//q + c'(u — w*) 1///3 + o((rt — w*) 1 /' 3 ) 


and the relation r e fj = dap/da, we deduce that the exponent 7 of (5.5) is given by 7 = 
1/(3 — 1/a, and thus a; = 07 is 


cu = 


2 +2k 
11 - k 


(5.9) 


This is a rational number which differs from (5.3). This implies that despite its appealing 
features, the Argyres-Douglas ansatz (5.1) cannot be exact. The problem of finding the general 
form of the /3 function for these non-local field theories thus remains open. Maybe one can guess 
the result by using the exact solutions for r e ff presented in this paper around the superconformal 
points and use the arguments above to relate r e fj and r. 

A simple inspection of the formulae (5.7) and (5.9) reveals that the two exponents a and 
to are related by 

u = 2(a-l). (5.10) 


Such a “scaling” relation is possible because, due to N = 2 supersymmetry, the only free 
parameter is k and thus one expects only one independent critical exponent, ft is actually 
possible to understand (5.10) on general grounds, following the line of reasoning of [7]: if we 
denote by U the U(l) N = 2 vector superfield contained in the low energy effective action at a 
generic point on the Coulomb branch, a deviation of the coupling constant from its fixed point 
value corresponds to the irrelevant (since the fixed point is IR stable) operator 


5L eS = A~“ 


d A eu 2 , 


where the fermionic integration is performed over half of N = 2 superspace. As the dimension 
of L e ff must be 4 since the action is dimensionless, and the dimension of d 4 d is two, we must 
have —uj + djj 2 = 2 where djj 2 is the dimension of the operator U 2 . Because of superconformal 
invariance, we must have djj 2 = 2 du = 2 a because these scaling dimensions are directly 
proportional to the R charges of U and U 2 under the exact quantum U(l)p symmetry of the 
superconformal algebra. We thus recover the relation (5.10). Conversely, one can consider 
(5.10) as a non-trivial test that the theory has indeed N = 2 superconformal invariance. 
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APPENDIX A 


For reference, in this appendix we briefly discuss the positions of the singularities and their 
flows for the massive Nf = 1,2,3 theories. For simplicity, all non-zero bare masses are taken 
to be equal. 


Nf = 3 with mi = m 2 = m 3 = m 

For m/0 the flavour symmetry is SU(3), while for m — 0 it is spin( 6 ) ~ SU(4). For the 
present case, the discriminant of the cubic polynomial ( 2 . 1 ), ( 2 . 2 ) in x is 


A \ o A 3 

A = — — U + 77T -I— -m 

16 


x 


u + -Asm 


4 ) “ “ A3m3 


256 


A 3 m 2 


2048 


A 


(A.l) 


showing that there is a triple singularity at, transforming as the 3 of SU(3), and two singlet 
singularities a± at 


2 A 3 

at — m H- m 


a± = -—A 3 m + 
16 


Ai 

512 


± A 


1/2 


m + 


M 

64 ) 


3/2 


(A. 2 ) 


' A 2 Q 

For small m one has a + — ^ and a_ ~ — IA 3 m, and at m = 0, a_ and at coincide and 
one has a quadruple singularity at u — 0 corresponding to a massless dyon of magnetic charge 
one which transforms as the 4 of SU(4). The singularity at a+ is due to a massless dyon of 
magnetic charge two [ 2 ] (see Fig. 2 ). As m is increased, the singlet cr_ moves to the left and 
the triplet moves to the right. The singularities at and a+ meet at 



Ai 2 
at = a+ = 32 = 2m 


(A.3) 


which is the superconformal point. As m is increased further, as discussed in Section 2 , the 
quantum numbers at the singularities are changed and now are ( 1 , l)o and (— 1 , l)o at cr_ and 
0 + and a quark triplet at at. The latter disappears to infinity as m —> 00 , A 3 —> 0, m 3 A 3 = Aq 
fixed, while a± —> ±Aq, and indeed, (after rotating the quark cut to the right and shifting 
a,D —> a,D + a ) we are left with the pure gauge theory Nf = 0. 


Nf — 3 with mi = m 2 = 0 and m 3 = m 


A 2 

For this case, the roots of the cubic are at x = u and x = ± 

i_ \ 

Thus there is a singlet singularity at a s and two doublet singularities at a^\ 


Ai (Ai _ 1 4 rr b 

256 l 64 + 


n 1/2 


a s — m + 


Ai. 

256 


± , a 3 

ai = ±—m . 
a o 


(A.4) 


As m is increased from 0 , the quadruple singularity at the origin splits into two doublets of 
massless monopoles (having different values of s) while the singlet at a s starts to move to the 
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right. But moves faster to the right and meets a s at the superconformal point 





Ai_ 

128 


(A.5) 


As m is increased beyond a s remains larger than a^ and the singularity at a s is now due to 
a (singlet) massless quark, while the other two doublet states at crj~ and a^ now are a magnetic 
monopole and the dyon with n e = =pl both having vanishing s. As m —> oo, 777 .A 3 = A| hxed, 
the massless quark at a s disappears to infinity, while aj —> ±-^, and one is left with the 
massless Nf = 2 theory. 

Nf = 2 with mi = m 2 = m 

This case was discussed in great detail in the main body of this paper, see Section 2. 


Nf = 2 with m\ = 0 and m 2 = m 

In this case the flavour symmetry is only U(l) x U(l) for m 7^ 0. The discriminant of the 


cubic is 


A A 2 

A = — 
16 


a 9 q A 9 9 9 a 9 Ao 27 t a 

u - m^u 6 - —u H-A u + ——-Afm 4 

32 64 2 4096 256 2 


(A.6) 


For m — 0 we have two doublets of massless monopoles and dyons in one and the other spinor 

A 2 

representation of spin(4) at u — ±^-. For small non-zero mass m each of these doublets splits 
into two singlets at 


= — — ± — m + 0(m 2 ) , crj: = — ± i—m + 0(m 2 ) . 

8 2 8 2 


(A.7) 


-*-y 

G 


-X- 


i 

x 


+ 



G 


+ 


Fig. 19: The flow of the singularities on the complex u-plane for Nf = 2, mi = 0 
and small m 2 = m 


Two singularities move off the real axis into the complex plane, see Fig. 19. As m is 
increased at moves to the right on the real axis, much faster than the other singularities. 
It passes between a _£ and oqi = (a [£)*. As m —> 00 , at —> 00 , while aZ and crj: arrange 
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themselves on a circle at the roots of u 3 + j^A^m 2 = 0. Since mA 2 = A 3 these are precisely 
the singularities of the massless Nf = 1 theory. In the present case, while flowing from 
m — 0 to m — oo, we never encounter a superconformal point where singularities coincide. 
Nevertheless, the quantum numbers of the massless states change by the monodromy matrix 
of another singularity as explained in Section 2. 

The Nf = 1 theory 
The discriminant is 

A = — —u 3 + m 2 u 2 + -A 3 mu — A 3 m 3 ——A® . (A. 8 ) 

16 [ 8 1 1 256 X J y J 

The positions of the singularities are given by the three roots of this cubic polynomial in u: 

uo = a + + a- + , a± = e ± 2 * 7 r// 3 o;_|_ + e T 2 wr ' / 3 a;_ + . (A.9) 

o o 

The massless states are all singlets. The a± are given by 



Fig. 20: The flow of the singularities for Nf — 1. The left figure shows the 
positions of the singularities for small mass m, the middle figure shows them at the 
superconformal point where <r + and <r_ coincide, while the right figure shows them 
for large m. 

While the choice of phase of the square root does not really matter since one can always 
exchange the roles of a .|_ and a_, the phase of the cubic root has to be determined carefully, 
so that the three singularities o are either all real or one real and the other two complex 
conjugate, and such that the singularities vary continuously as m is varied from 0 to oo. In 
particular, for 0 < m 3 < §|A 3 the curly bracket in (A. 10) is real and negative, and we choose 
a± = — i|{.. •}| 1 ^ 3 - This ensures that ctq is real and negative, while < 7 + and <r_ are complex 
conjugate. For m — 0 one simply has ctq = — uq , cr± = e ±wr / 3 wo with uq = 3A 2 /(4 • 2 2 / 3 ). As 
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m is increased, <to moves to the left and a± approach the real axis, see Fig. 20. For small m 
one has 

ao ~ — mo(1 + 2 n) , (T± — uq ^e ±i7r//3 + 2pe =R7r//3 j (A.11) 

where fi = 2 ■ 2 1 / 3 m/(3Ai) <C 1. When the square root in (A.10) vanishes, a + = cc_, so that 
<y+ and <T— meet on the real axis. This is a superconformal point: 


m = 


3 Ai 



3Al 

~4~ 


(A.12) 


As m is increased further, the curly bracket in (A. 10) is now complex. We choose [.. ,] 1//2 = 
—i |[.. ■]| i/2 . (The opposite choice only exchanges <r + and <r_.) The phase of the cubic root is 
chosen so that a+ = a*_ with ^ < arga + < tt and —it < argct_ < —yf- With this choice, 
the singularities will flow continuously: <to keeps moving to the left, while <j±, now both real, 
move to the right, but <t+ moves faster, see Fig. 20. As m —> oo, one has 


do — —(Afm) 1 / 2 , (T_ ~ (Afm) 1 / 2 , <r + ~ m 2 (A. 13) 


up to terms 0(A\/m). Since A 3 m = A|, cro and <r_ flow to the singularities of the pure gauge 
theory, while cr + disappears to infinity. 


APPENDIX B 


(j) 

In this appendix we will show how to express the integrals /• in terms of the standard 
elliptic integrals K(k), E(k), Ui(v,k), and also give their expressions in terms of the Weier- 
straB p function, as well as give certain relations between elliptic integrals. Basic references 
are [14, 19, 20]. Our conventions are those of [14], Unless otherwise stated, we will assume 
that the roots ej of the cubic r] 2 = 4£ 3 — £ — (ft are all different. They obey ]>T = 0, and 

V 2 = 4 n,(£-o)- 

The K , E and III only depend on the square of k which was defined in terms of the ej as 
k 2 = t 2 ~t 3 • We use the notation* of [14]: 

ei—e3 i J 


K(k) 
E(k) 


III (u, k) 


dx 


[(1 - x 2 )(l - k 2 x 2 )] 1/2 


dx 


1 — k 2 x 2 \ 1/2 


1 — x 2 


dx 


[(1 — x 2 )(l — k 2 x 2 )] 2 ^ 2 (1 + ux 2 ) 


(B.l) 


-k Note that in Mathematica which we used extensively for the numerical determination of the curves 
of marginal stability, one denotes A'(fc) = EllipticK[fc 2 ], E(k) = EllipticE[fc 2 ] and Hi(v, k) = 

EllipticPi[ — v, k 2 ]. 
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We will now show that the integrals over the cycle 71 are 


where 


/«=2 


V (ex-es) 1 / 2 


e3 

e2 


/<'> = 2 1 ^ 


V (ei-e 3 ) 1 /2 


K(k) 


[eiK(k) + (e 3 - ei )E(k)] 


e3 


62 


d£ 


4 1} = 2 

J v(£-c) 


e3 


(ei - e 3 ) 3 / 2 
4 k' 


1 — c + k' 
1 


AT(&) 


1 + k! (1 - c ) 2 - k' 2 


ni u(c), 


1 - k' 
1 + k' 


k' 2 = 1 - k 2 


c = 


C - e 3 
ei - e 3 


u(c = 


1 -c+k'\ 2 f 1 — k' x 2 


1 — c — k'j \1 + k' 


(B.2) 


(B.3) 


( 2 ) 

The corresponding integrals /> ' over the cycles 72 are obtained from equations (B.2) and 
(B.3) by exchanging everywhere in these equations the roots e\ and e 3 . In particular, this 
exchanges k and k', and results in the exchange of c and 1 — c, so that u(c) gets replaced by 

_ ( c+k \ 2 ( 1 — k \ 2 
\c—k) Vl +k) ' 

To convert the Ij into the K , E, Ill one needs to transform the cubic curve into a quartic 
curve. Let’s demonstrate this for ij 1 "*: 


62 62 k 2 



f ^ 1 . 

( d(' 

1 V J 

[(£ - e i)(£ - e 2)(£ - e ^)} 1 / 2 (ei - e 3) 1 / 2 J 

[(#- i)A- fc 2 )d 1/2 


(B.4) 

where we changed variables , £ = (£ — e 3 )/(ei — e 3 ), and introduced k 2 as above. The 
transformation to a quartic curve is achieved by the further change of variables 

| = 1 + k' + —-— (B.5) 

^ ~ W 


where the complementary modulus k' is given by (B.3). (The choice of sign for k' does not 
matter.) A final rescaling x = 2k! C yields 



1 2 f 

(ei - e 3 )i / 2 rTA 7 J 

-1 


dx 


(MV 

V l+k') 


x* 


(1 


x “ 


1/2 


_ 1 _Li-fiA) 

(ei-e 3 ) 1/2 l + *' \l+k'J 

(B. 6 ) 
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Using one of the standard relations between elliptic integrals with different moduli [14] 


2 K f -—= K(k) , (l + k')E(l—% ) =E(k) + k'K(k) , (B.7) 


1 + k! V 1 + k' 


1 + k f 


one obtains the hrst equation (B.2).* 


(1) 


Going through exactly the same steps for the integral A; ; leads to 


/(B _ 


(ei — 03) 1 / 2 1 + k! 


dx- 



/ \ 2 \ -1 


(ei - e 3 ) 

1 + * 2fc'(i (J;*) x 2 ) 

+ e 3 


-1 


i-(hU 2 ^j(i-^ 2 ) 

l-k! 


(ei — e^) 1 / 2 1 + k' 


((ei - e 3 )(l + k') + e 3 ) A' 


1/2 


1 + k' 


>(-(££) • 


1-J/v 1-c 

1 + fc' 


(B. 8 ) 


Using the relation (1 — A; 2 )Uli(—/c 2 , k) = A(/c) with k = (1 — A; / )/(l + A: / ), as well as the relations 
(B.7), one arrives at the second equation (B.2). 

Finally for I^\c), going through the same steps leads to 

1 



1 2 r 

(ei - e 3) 3 / 2 1 + k' J 




(B.9) 


1 

1 + k! — c 



2 k' 

1 + k! — c 


1 


1— fc' y, l — k' — C 

1 +k" L l+k'-c 


where c = (c — e 3 )/(ei — 63 ). The last term in this expression is of the type a J_ b = and 

can be replaced by a z i; z _ b 2 ■ It is then clear, using again (B.7) that one gets the third relation 
(B.2). . 

The integrals I%\c) can be simplified if c is one of the roots e^. We will always consider 
4 1 } (c). Everything can be translated for I^\c) if we permute everywhere k and k' as well as 


-k Note that it is quite non-trivial to keep track of the correct overall sign. For example, to get the last 
equality of eq. (B.6) one needs to define carefully where the cuts of the differents square roots lie. For 
real masses m, and real A, all one can get is a sign ambiguity that may depend on sign(3?nit). This 
ambiguity is the same for all three integrals Ii, I 2 and I 3 , so that it is most easily fixed by comparing the 
resulting a and an with the required asymptotics (2.25). 
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e 3 and e\. Let first c — e Then c = 1 and v W = - (w) = ~~ k2 (see (B.3)). One has [14] 


ni (~~k 2 ,~k) = 


1 i _i_ y 

2 “ - E(k) = -4n~ (E(k) + k'K(k)) . 


1-k 2 


4k' 


Hence 


4%i) = 


E(k) 


(ei - e 3 ) 3 / 2 A; 2 - 1 

r(l)/_ \ _ o d rU) T Toin rr Z-2 d TS ( 1 „\ _ 1 


(B.10) 


(B.ll) 


This can also be directly obtained since ( e \) = 2^/j . Using k 2 -^ I K{k) = —(/c) + 

this gives again (B.ll). Obviously, this constitutes a consistency check for the third 
equation (B.2). Let now c = e 2 - Then c = /r and z/(c) = —1. Hi is singular for v = —1, but 
I^\e 2 ) can be obtained in exactly the same way since I^\e 2 ) = 2 ^/] ^ provided we keep 
the integration cycle 71 fixed and away from e 2 (and e 3 ). We then immediately get 


Ute) = 


gw 

(d -e 3 )3/2 A ; 2 Vl - k 2 


- K{k) 


(B.12) 


Finally, if c = e 3 , one has c = 0 and also z/(c) = —1, and again Hi is singular, but we can still 
use I^\es) = 2 ^/j 1 ^ which readily gives 


Uw) = 


(ei - e 3 ) 3 / 2 k 2 


73 WW - SW) • 


(B.13) 


Of course, using the first two equations (B.2), one can reexpress K(k ) and E{k ) in terms of 
/j 1 ^ and and hence we have shown how to express I^\ej) in terms of and I^\ 

(j) 

We also need to Taylor expand ig (c + 5c), hence we want to compute the derivative 
^ lip (c). For any cycle 7 j we have 

d£ 


dc 3 {C) ~ 


v(£ - c ) 2 


73 


(B.14) 


Now observe that 


0= (f d( 


V 


d ^ - c 


= $ d£ 


12g ~ 92 
2 ??(£ - c) 


d£ 




(B.15) 


This is a sum of integrals containing /i ,/2 and J 3 , as well as the integral (B.14). Solving for 
the latter gives 


(4c 3 - C0 2 - 0 3 ) 


d£ 




= —2c 1 1 + 2/2 


i(12c 2 -0 2 )/ 3 (c) . (B.16) 


Inserting now the appropriate values of 02 and g 3 for the Nf = 1 curve and taking c = — u/3 
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(B.17) 


yields the desired relation for the derivative /3(c) = g^/3(c): 


32 


+ h- 


A word of caution is in order: in eq. (B.2) we have replaced the integrals over the cycle 71 
by twice the integrals from e3 to ■ Depending on the detailed form of the cycle, these two 
definitions may differ for I3, which has a pole at £ = c, by terms 27rires^-. So we allow the 
freedom to add such terms “by hand”, this being equivalent to changing the definition of the 
cycle 71 with respect to the position of the pole. A related point is that for c = e2 or c = ez, 
one has z/(c) = —1 and III (—1, k ) diverges. Nevertheless, we can remove the divergence by an 
appropriate choice of integration contour! With this choice fli (—1, k ) should be understood as 
K(k) — E(k)/( 1 — k 2 ). Then, eqs. (B.12) and (B.13) also follow in a straightforward manner, 
using (B.7), from the third equation (B.2). Similar remarks apply to the integrals /■ ; over 
the cycles 72. 



Fig. 21: The elliptic curve in the z-plane where it is simply a parallelogram with 
sides cai and u>2- The cycles 71 and 72 are straight lines parallel to the sides of the 
parallelogram. 


To gain a somewhat better control over the integration cycles 71 and 72, in particular with 
respect to the positions of the poles, it is sometimes advantageous to introduce the uniformizing 
variable z via the doubly periodic Weierstrafi p function [19] as £ = p(z), 77 = p'(z). Then the 
integrals l[‘^ are simply 



7 i 


f p\z)dz 

J 0(z) 


LUi 

J dz = Ui 
0 


(B.18) 


where and u>2 are the two periods of p and are such that e\ = p (77-), = p ( 1 ^ 1 + LJ2 ), 


f This is so because v —> 1 typically occurs under the RG flow as some bare mass rrij —> 00 and the 
divergent part of IIi precisely is some integer multiple of 777=. 
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e 3 = P (if)- The cycle 71 is mapped in the z-plane to a straight line from 0 to 07 (or any 
shifted copy of it) while 72 is mapped to the straight line from 0 to CJ 2 , see Fig. 21. 

Obviously one has r = Given the ei,e 2 ,e 3 , one can obtain ccq and u> 2 , hence /j 1 ' 
( 2 ) 

and I\ using the inverse of the p function. The latter is conveniently expressed through the 
inverse of the Jacobi elliptic function sn: 




(ei - e^) 1 / 2 
r (0 


sn 


-1 


ei - e 3 \ 
£ - e 3 ) 


1/2 


,k 


(B.19) 


For a numerical computation of the c <7 — l[ , eq. (B.2) is of course simplest. Next one has 
(do not confuse the rji with rj) 


Ai) _ 
1 2 — 




= -2 m 


(B. 20 ) 


7* 


where again numerical computation of the ry is easiest done via (B.2). Finally, one has 

1 


r(*) _ 
— 


- c) p'{z C ) 


[2ca*C(^c) - 4^27 


(B.21) 


7 i 


where z c is given by (B.19) with £ = c, and 


: / . 2m 7 t 9\ (nz 

C(z) = —z + —^( — 

UJ\ 0J\ \ Oil 


(B.22) 


APPENDIX C 


In this appendix we give the differentials A for the different massive theories as well as the 
decomposition of the period integrals in terms of the three basic elliptic integrals I\, I 2 and 
J 3 . We also check some RG flows at the level of these integrals. 


One flavour: Nj — 1 

A one-form A satisfying (2.4) obviously is 


A = 


y/2 y dx 
47T x 2 


\/2 

47T 


dx 


dx 


— \ 3x — 2u + 


m A dx 


4 x J 2y 


(C.l) 


The first term in the bracket is an exact form and vanishes upon integration over a cycle. 
Converting to Weierstrafi normal form (2.21) by 7 = 2 y, £ = x — ^ we get 




A? 


u I\ - 3 1 2 - —mh -o 


u 


(C.2) 


We see from (C.l) that A has poles at (x — 0, y — ^iAf/8) with residues in agreement 

with eq. ( 2 . 10 ). 
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Two flavours: Nf — 2 
We have [ 2 ] 


A 


y /2 y dx 

47r r 2 _ M 
x 64 


\/2 dx 
47 T y 


x — u — 


\/2dx4y 2 f 1 1 \ 

4 tt y A| x + fj 

A| (mi — m2) 2 i A| (mi + m2) 2 

16 r _ M. + 16 “ , M 

x g x W 8 


(C. 3 ) 


Converting to Weierstrafi normal form of the cubic, again by 77 = 2 y, £ = x — | we arrive at 


A = 


V2 

47 r 


4 T t A9 .o t /AoU 

-// /, -2/2- y(mi -m 2 y h f Y - 3 


An , \2 T ( ^2 ^ 

T (ro i + m2) *(“-3 


(CU) 

One sees from (C. 3 ) that A has poles at (x = ^,7/ = iiA^mi — 7712 )/ 8 ) with residues 
and at ^ = -¥>y = TiA 2 (mi + m 2 )/ 8 ) with residues ^ In 

particular, for mi = m2 = m, there are only the poles at (x = — -^,7/ = TiF^m/ 4 ) with 
residues ±5 


Three flavours: Nf = 3 

This case is more complicated since the y 2 is no longer linear in u, see eq. ( 2 . 1 ), ( 2 . 2 ). 
There exist expressions for A in the literature using a quartic curve [ 21 ] instead. Proceeding 
along the same lines as in [21] we can obtain A also for the cubic curve (2.1), (2.2): we write 


y 2 = G(x) — F 2 {x) 


F{x) 

G(x) 


x 
4 a 


a yx — u 

4 b 2 x 2 6 4 
2 a 4 a 


x“ . 6“A 


2 a 2 a J 


ex 


(C. 5 ) 


where we have set a = b 2 = a nr 2 , c = |A3 ]Q m* and d — a ^ 2 i<j mfm 2 . One can then 
check that the following differential indeed satisfies ( 2 . 4 ) 


A 


V^xdx /F(x)G'(x) , 

877 y/ay \ 2 G(x) 


(C.6) 


In general, to express J A in terms of the three elliptic integrals Ij one needs to decompose 
the quartic polynomial G(x) into linear factors. While this can be done in general, it is 
cumbersome and the result not very illuminating. Here we will restrict ourselves to the simpler 
case m\ = m2 = 0 , m3 = m. Then, since c = d = 0 : G(x) = ^(x 2 — b 2 ) 2 and ^ 
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where now b = Asm/8. It is then easy to see that 


A - 

8tt y 


x 


n b 2 + bu b 2 — bu 
2 u +-— +-r- 


x + b 


x 


(C.7) 


so that after introducing rj = 2 y and £ = x — ^ we get 


A = 


V 2 

47T 


Aq 


-u — 


192 J 




64 

A3 m 
64 


( 8 m + A 3 m )/3 


( 8 m — A 3 m )/ 3 


Asm 


Asm 


A 2 

A 3 

192 


A| 

192 


(C.8) 


We see from (C.7) that A has poles at {x = A 3 m/8, y = =p(A 3 m/8 — m)A 3 /8) and at 
(x = —Asm/8 , y = ±i(—A 3 m/8 — m)A 3 /8) with residues Note that the integrals J 3 

cancel in the m —> 0 limit. 


The pure gauge theory: iVy = 0 

Finally let us note that for the pure gauge theory in the conventions of [2] we have 


/ 


A 


V2 

47T 


A 4 

M/l - 3/2 - 



(C.9) 


Since — | is one of the roots e* one can use relations (B.ll) to (B.13) to reexpress J 3 (—|) as 
a combination of /1 and I 2 ■ This must be so for Nj — 0 since A has no poles and thus its 
integral must be expressible through I\ and I 2 only. 


RG flows of the integrals 

One can check the different RG flows at the level of the period integrals. For example, 
starting with the Nf = 1 periods (C.2) and letting m —> 00, Ai —»• 0 while keeping mA'\ = Aq 
fixed, we immediately find that (C.2) flows to (C.9). The flow from Nf = 2, eq. (C.4), to 
Nf = 1, eq. (C.2), as m 2 —> 00 , A 2 —>► 0, m 2 A 2 = A' 4 fixed, is less trivial: starting with (C.4) 
we have 


Nf= 2 


V2 

47T 


4 A? 

-m/i - 2/ 2 + — -(m 2 - 2m 1 )Is 
o o 

A 3 

- + 2mi)Is 


u A\ 
3 87712 


M 

’3 


AL 

87712 


(C.10) 


Here, the integrals J 7 ; on the r.h.s. are meant to be those of Nf = 1. Taylor expanding 
J 3 (c + Ac) = J 3 (c) + J 3 (c)Ac + ..., and using the relation (B.17) for the derivative of J 3 (valid 
for IVj = 1), the r.h.s of eq. (C.10) becomes exactly the r.h.s. of eq. (C.2), up to terms that 
vanish as m 2 —>► 00. Similarly one can check that as m —> 00, the = 3 periods (C.8) flow 
to the Nf = 2 periods (C.4) with 7771 = m 2 = 0. This requires to reexpress Is (c) for c a root 
e,; in terms of Ji and I 2 through the formulae given in appendix B. 
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APPENDIX D 


In this appendix we perform some checks on our equations (2.27) that express a and ajj 
for Nf = 2, mi = m2 = m in terms of the 1 / and hence of the complete elliptic integrals K , 
E and III. In particular, we will show that one indeed recovers the correct expressions of the 
massless Nf — 2 theory [4] as m -> 0, and that one gets the appropriate expressions of Nf = 0 
in the m —> 00 limit. 


First we examine the limit m 
has 


0. Then the extra term 


m 

s /2 


111 


(2.27) disappears and one 



7 i 




V2 r df (t« - 2g) 

4vr J 77 

7 i 


s /2 / d x(u — x) 

An J y 

7 i 



7 * 


dxs/x — u 



(D.l) 


2 A 2 a 2 

Since ei —» 62 —> — § + e3 —» — | ^ (or in terms of the variable x the roots are 

A 2 

xi = u, X 2 = -/), these are precisely the integral expressions for ajj and a of the massless 
Nf = 2 theory (which equal ^ times the ajj and a of the Nf = 0 theory) [1, 4], 

Now let m / 0. It is then straightforward to explicitly check that at the rightmost 
singularity, u = m 2 + we have the following: if m < // one has e\ = e2, k! — 0 and 

thus it follows that ao(u) — 0; if m > however, one has 62 = e3 so that k — 0 and it 
then follows that a(u) = -^=. Hence we see that the massless BPS state at this singularity 

u = m 2 + is a magnetic monopole (n e ,n m ) s = (0, l)o if m < while for m > it is 
a quark (n e ,n m ) s = (l,0)i. This perfectly agrees with the discussion of Section 2 and thus 
justifies the choice of adding the -^= term in the first equation (2.27). 

Let us now examine the RG flow to the Nf = 0 theory as m —> 00, A2 —>► 0, 777.A2 = Aq 
hxed. The flow of the roots e* of eq. (A.9) is 


ei 


62 


e 3 



(D.2) 


where eis the standard labeling of the roots for Nf = 0. We see that the flow exchanges 
the labelling of e2 and 63. Let us first consider the flow of a(u) as given by the first eq. (2.27). 
Since the Nf = 2 curve flows to the Nf = 0 curve, the integrands flow appropriately. Clearly, 

the cycle encircling e2 and e3 flows to the cycle r y[ Nf ~°\ Hence —> /^|jv /= o and 

I^ —■> 1 ^ |_/V/=o- integral /g 1 ^ is more subtle since it is the one that involves the pole at 
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A 2 (0) A 2 (0) A 2 

c = — I —jf = e 3 —jp As A2 —0, this pole approaches eg . Hence, for non-zero ^ the 

pole at £ = c is outside the integration contour. As A2 —> 0, the pole crosses the contour and 

/g 1 ) picks up a contribution from the residue which is 51^ = which precisely cancels the 

additional term -^= in eq. (2.27). Hence 


am 


V2 

47T 


-ul[ 1] -21^ - ^0/. (1) 
3 1 2 2 


fe (0) ' 
3 l e 3 > 


(D.3) 


Nf =0 


Using the relation 


Nf =0 


r(l) 


+ ¥/■ 


(1) 


V J 2 ^ 3 1 


Af /=0 


from appendix B, we see that the 


r.h.s. of eq. (D.3) indeed coincides with the corresponding integral for Nj = 0, and hence 
= cT°)(w). Next, let us discuss the flow of od(w). The cycle 72 encircles 


am 


am 


Nf=0 


e\ and e2 which flow to e[ (l) and 


The corresponding integral thus is the sum of the 


integral over a cycle 7^ around e^’ and e^’ and one over a cycle 7) u ' around e^’ and eg 

72 —> 72 0) + e 7i°^ where e = sign(Am'u). Again, the pole crossing the cycle yj 0 '* gives a term 
—e^=. As a result we have 


,( 0 ) 


To) 


(0) 


To) 


TO). 


a(u) —■» a,(°\u) 

a D (u) + e —= —> a,p\u) + ea^(«) . 


(D.4) 


This motivates us to define d]j(u) = ajj(u) — e ^a(u) — -^=J which is such that under the RG 
flow as m —» 00 one has a(u) —> gT°)(m) and do{u) —> a^\u). 
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